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Abstract 

We derive the angular correlation functions between an arbitary 
number of partons inside a quark or gluon jet emerging from a high en- 
ergy hard collision. As an application results for the correlation in the 
relative angle between two partons as well as the multiplicity moments 
of any order for partons in a sidewise angular region are calculated. 
At asymptotic energies the distributions of rescaled angular variables 
approach a scaling limit of a new type with two redundant quantities. 
All observables reveal the power behaviour characteristic to the fully 
developed selfsimilar cascade in appropriate regions of phase space. 
We present analytical results from the double logarithmic approxima- 
tion as well as Monte Carlo results. 
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1 Introduction 



Originally Perturbative QCD has been applied to the calculation of quantities 
with the hadronic final states largely averaged over, either total cross sec- 
tions or jet production in hard collision processes. Such calculations turned 
out to be very successful in comparison with experiment. On the other hand 
there is continued effort to apply PQCD also to the more detailed descrip- 
tion of the final states [|l], |^. Such studies have been performed in Double 
Log Approximation (DLA) summing the infrared and collinearly divergent 
terms and taking also into account the soft gluon interference effects by "an- 
gular ordering" 0, • A prominent example is the prediction of the energy 
dependence of total multiplicity in a jet 0. Another noteworthy result is 
the derivation of the "drag" or "string" effect suggested earlier on phe- 
nomenological grounds 0. Also the one-particle momentum spectrum is 
successfully predicted indicating that the QCD evolution can be meaning- 
fully continued to low momentum scales of the order of hadron masses which 
led to the notion of "Local Parton Hadron Duality" (LPHD) p|. Recently 
the two-particle momentum correlations have been studied ||^ and the shape 
although not the normalization had been experimentally verified |10|. The 



aim of such studies is not primarly a further test of PQCD at a fundamental 
level, rather one wants to find out the limiting scale for the application of 
PQCD, and thereby to learn about the onset of nonperturbative confinement 
forces. 

This paper summarizes our study of further details of the final states in 
QCD (for selected results see also [^-[|T^). We give a systematic derivation 
of a variety of multiparton correlations starting from the master equation for 
the generating functional in DLA |T^. This approach allows to attack for the 
first time the computation of observables for multiparticle final states and 
to find some simple asymptotic formulae. All the observables considered ap- 
proach a universal scaling limit at high energies described by a single function 
of a scaling variable. In this way quite different observables become closely 
related. The validity of our analytic formulae is checked by Monte Carlo 
calculations (using the program HERWIG |jl6|)- The experimental confirma- 
tion of such relations would support the idea that the underlying structure 
of multiparticle production is the parton cascade with little disturbance by 
confinement effects. 

Our results also allow to directly adress the question of the fractal nature 
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of the QCD cascade [|T^, Indeed we find the kinematic regime where 
our observables depend hke a power on the resolution scale. This behaviour 
is found sufficiently far away from the cutoff scale. The simplest exam- 



ple is the well known result for the total multiplicity n ^ exp(/5yP0/A) = 
[PQ / A)"-^^^^ in the cone with angle 6 of a jet with momentum P. The energy 
dependence is determined by the anomalous dimension a = 70 = ^JGoJti. 
For the higher order correlations we find power laws of similar type in the ap- 
propriate phase space region. The power behaviour of multiplicity moments 
suggested as consequence of an underlying scale invariant dynamics ("In- 
termittency" |19|) has been looked for in many collision processes pO). We 



specify the variables and the kinematic regime appropriate for such studies 
and calculate explicitly the intermittency exponents for the QCD cascade. 

The problem of angular correlations with application to azimuthal corre- 
lations has also been treated in detail in [21| where the same scaling property 
at high energy is found as in case of our variables. Results on multiplicity 
moments have been obtained in parallel works by two other groups ^ 
applying the method of which is based on the same principles but dif- 
fers in details. Whereas also gives some corrections due to the MLLA 
and [^] derives further phenomenological applications, we first construct the 
fully differential n-particle angular correlations and then find the universal 
properties of various specific projections. 

The paper is organized as follows. In the next section we derive the inte- 
gral equations for multiparticle distributions of arbitrary order and discuss 
our approximations. In Section 3 we give results on single particle observ- 
ables partly known already. Section 4 explains our methods to solve the 
integral equations for multiparticle angular correlations. These are applied 
in Section 5 to the derivation of the two parton angular correlation func- 
tion, the n-particle cumulant correlation functions and in Section 6 to the 
n-particle multiplicity moments in different angular regions. In Section 7 we 
compare our analytic results with Monte Carlo calculations. Conclusions are 
presented in the last Section. The Appendices contain various cross-checks 
and details of the calculations. 
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2 Integral equations for multiparton correla- 
tions 



2.1 Generating function 

A compact description of the multiparticle distributions is provided by the 



generating functional [24 



Z{u} = J2 f dkl... f dklu{ki) . . . u{kn)p{ki, ...,kn) (1) 

where p{ki, . . . , kn) is the probabihty density for exclusive production of par- 
ticles with 3-momenta ki . . . kn and u{ki) can be thought of as profile or 
acceptance functions. Z{u} contains the complete information about the 
multiparticle production. For example ^6"'/6u{ki) . . . 6u{kn)Z{u}\u=o = 
p{ki,...,kn) reproduces the probability densities, while Taylor expansion 
around u{k) = 1 generates the inclusive densities 

p(")(A;i, = 6^Z{u}/6u{k,)..Mkn) \u=i • (2) 

Setting u{k) = v for all k gives the generating function of the distribution 
of the total multiplicity. Let us also introduce the connected correlation 
functions (also refered to as cumulants) which are obtained from 

r(")(fci, ...,kn)= 6^\nZ{u}/6u{k^) . . . 6u{kn) \u=i • (3) 

They are related to the inclusive densities by 

pH ^ + 4^),. (4) 

where (ip^o^ is build from products of correlations of lower order than n, for 
example 

<1(1,2) = P^'\1)P^'\2) (5) 
41(1>2,3) = np«W+p«(l)r(2)(2,3) + ci/c/. (6) 

4 = 1 

The generating functional approach to QCD jets has been first developped 
in Refs.0, In the double log approximation (DLA) the master equation 
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for the generating functional of multiparton densities for an initial parton a 
(either a gluon g or a quark q) reads |jl5 



Zp,a{u} = exp( [ Mp,a{K)[u{K)ZKAu} - l]d^K] 
\jrp{K) j 



(7) 



The subscript P denotes collectively the momentum vector of the parent 
parton and the half opening angle [P = {P, 0}) of a jet it generates; Tp{K) 
stands for the phase space of the intermediate parent K {rp{K) = {K : 
K < P, Qkp < ©5 KQkp > Qo} where Qq is a transverse momentum cutoff 
parameter). Aipa{K) is the probability for bremsstrahlung of a single gluon 
off a primary parton a which reads for small angle Qpx 

Mp,a{K)d'K = Caa\KQpK) — -r^—^, (8) 

^PK ^71" 

where Cg = 1, Cg = Cp/N^ = 4/9 and = = Gas/n is the QCD anoma- 
lous dimension controlling the multiplicity evolution. We treat either as 
constant or as running with a?{pT) = Z^^/ ln(pT/A) = /3^/[ln(pT/Qo) + 
with the QCD scale A, A = ln(go/A) and /?2 = 12(^N^ - fAT/)"^ (we use 
Nf = 5, so = ||). In this approximation the recoil of the radiating parton 
is neglected, i.e. its momentum remains unaltered. Also in (|^) the pair cre- 
ation of quarks inside a jet is neglected (the production of gluons dominates 
because of an infrared singularity and the larger color factor). 

Densities defined by Eqs. (@J^) describe only the radiated gluons excluding 
the primary parton. For example, the total multiplicity n = J p^^^d^K does 
not contain the beam (leading) parton which also contributes to the final 
state. To account for the leading particle one can define another generating 
function |^ 

Zp,aM = u{P)Zp,,{u}. (9) 
The master equation for Zp^a follows readily from Eq.(0) 

Zp,a{u} = <P) exp ( / Mp,a{K) \ZK,g{u} - ll d^K] , (10) 

\JTp(K) J 

For many problems, like production of particles into sidewise detectors this 
distinction is irrelevant. For fully inclusive observables the leading particle 
should be included; in cases studied by us with an angular average its effect 
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however decreases with increasing energy. We shall therefore mainly use the 
first formulation, pointing when necessary, to the differences between the 
two. 



2.2 Inclusive densities 

By differentiations (|],|D of the generating functional (0) we get the integral 
equations 

rPih, .., kn) = d^:L{K .., kn) + / (fKMp,a{K)p''g),{kr, .., k^) (11) 

where the "nested" term 

dplestih, k„) = Mp,a{ki)p^^ki\k2, .., kn) + cyd.perm. (12) 

refers to the direct emission of one parton followed by the "decay" into the 
remaining n — 1 ones. Here and in the following we drop the index a of 
correlation functions if the equation holds for both cases a = q and a = g. 
Using (^ we get integral equations for densities, explicitly for n = 1,2 

pfi^k) = Mp,aik)+ J^Mp,a{K)p^^^l{k)d^K, (13) 

P^p\ki,k2) = P^\k,)p'^^\k,)+Mp,a{kr)p';^l{k,)+Mp,a{k2)p^^l{k,) 

+ j^Mp,a{K)p%{ki,k2)d'K, (14) 

and for arbitrary order (see Fig.l) 

pf{k^,...,kn)=d!-p\ki,...,kn)+ j^Mp,a{K)p^;;},{kr,...,kn)d'K. (15) 
with = d^^Xd + dnesf 

Contributions containing Born terms M.p are non-leading in the high 
energy limit so they can be neglected at asymptotically high energies. For 
example, for n = 2 the "nested" term M.p{ki)p^ll{k2) , which corresponds 
to chain emission P — ^ fcl — > fc2, is nonleading at large P as compared to 
the product term pp\ki)pp\k2). Here, and in the following, the subscript 
r denotes generically all boundaries of the phase space integration. 
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Cumulant correlations containing leading particles are easily derived from 
Eqs.(BI) 

n 

ffik,, .., K) = {~ir-\n - 1)! n 5{P - h) + rPih, .., K) (16) 

1=1 

The corresponding densities are found from Eq. (P to be linearely related 
to the densities without leading particles, for example 

~p?{k) = 6iP-k)+p?ik), (17) 
pP{h,k2) = 6{P-k,)p^^\k^) + 6{P-k^)p^^\h)+pf{k,M){'^^) 

These quantities obey similar integral equations as in (p!3| , p!^ ). 

The above equations all hold for either quark or gluon jets. It should 
be noted, however, that in the jet we assumed production of gluons only. 
Therefore in the integral equations the correlation functions under the inte- 
gral always refer to gluon initiated intermediate jets whereas the quark degree 
of freedom only enters the emission from the primary parton described by 
P^a in our formulae. 

The difference between quark and gluon jets is seen most easily for the 
cumulant correlations which correspond to one primary gluon emission and 
are therefore proportional to Ca- From Eqs. (|Tl],|T2D one finds 

r£ = c.rg (19) 

From this result the more complicated equations for p^p\ or the correlations 
with leading particles can be derived. 

Finally we introduce our notation for inclusive densities. Unless stated 
otherwise a distribution p{x) is differential in its argument x, i.e. p{x) = 
dn/dx. For example, the distributions in spherical and polar angles are 
related by p(fi) = p{d)/{2'Kd). 

2.3 Connected correlation functions 

In the solution of (|T5[) a careful treatment of the boundaries F is required. 
These depend on the multiparticle kinematics and therefore also on the struc- 
ture of the direct terms d^p\ In particular, if one solves Eq.(|T5|) by iteration 
one finds in general different bounds for the first iteration (integral / d^Kd^^^) 
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and for all higher iterations. For illustration we consider in more detail the 
bounds for the angular integrals in the left column ("precise bounds") of 
Fig. 2. The simplest case is Eq. (p!3|) for p^p\^li). In the lowest order there is 
the chain P ^ K ^ 1. The limits of integration come from the lower px 
cutoff {kiOxki > <5o) the small circle around particle 1 in Fig.2a) and the 
requirement of angular ordering, namely, the intermediate parton K can emit 
parton 1 only in a cone around K limited by P, i.e. 6;^^^ < QpK, which 
leads to the outer boundary in the figure. For two particles the integral 
bounds / d^Kd^^^pj.^^ involving the product term have, in the first iteration 
of Eq.(|l4D, to obey the same constraints as above twice which yields the 
contour in Fig. 2b (bound F"). If we iterate Eq.(n) once more according to 



the chain emission P ^ Ki K2 (1? 2) the angular ordering requires the 
same outer bound for Ki as for K2 before but now the lower limit in the Ki 
integral is given by the minimal circle enclosing the final particles, because 
of angular ordering Bi^jfci, QK2k2 < ^KiK2 (see Fig.2.c for the i^i-integration 
domain). Similarly, bounds in all higher iterations are sensitive only to the 
minimal virtuality of the last parent (bound F'). 

This distinction applies to an arbitrary number of partons in the final 
state: whereas the bound F" depends on the coordinates of all individual 
momenta, the bound F' depends only on those of the circle. It is then ad- 
vantageous to work with the connected correlation functions. They satisfy 
the following equation, c.f. (§,|11|) 

^p}i{^lj--}^n) = C^P,7iest(^l; ^n) + ^P,l,(^l; --J ^n) 

+ j^d'KMp,a{K)V^^},{h,..,K) (20) 

which has the same boundary F' for all iterations. On the other hand the 
inhomogeneous term is defined in terms of the different boundary F" 

Agi(fcl, .., K) = d'KMp,aiK)d^^kproAkl, .., K). (21) 

The contribution of the nested term can be neglected in the high energy 
limit. 

It follows from Eqs.(pO|,21) that the connected correlations contain one 



more iteration of the Born kernel than the corresponding densities. This has 
two consequences: a) the perturbative expansion of the cumulants starts at 
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one power of a? more, and b) the connected and disconnected correlations 
have different factorisation properties. This last feature will be exploited 
later in detail. 

2.4 Pole approximation 

The double logarithmic approximation consists of identifying regions of the 
phase space which give dominant contributions to the integrals involved. 
In these domains only simple poles from inner emissions dominate - all 
other slowly varying functions are replaced by their central values at sin- 
gular points. This "pole dominance" approximation is an important part of 
the DLA and was used by many authors for calculations of the total multi- 
phcities and momentum distributions 



2.4.1 Single parton density 

As a first example we solve Eq. (|13D for the single parton inclusive distribution 
in k. Writing the momentum integration explicitly for dn/dVtdk with k = \k\ 
(see Fig. 3), 

p«(^], k) = k'Mp{k) + ^ f ^^a\KQpK)pi\nKk, k). (22) 

The dominating singularities are those in Qxk, ^Kk ~ since QpK > ^Kk 
due to the angular ordering. At this point Qpx ~ ©Pfc = 'd- Azimuthal 
integration around the k direction is trivial and we get 

^KkPK\^Kk,k)] , 

(23) 

where the factor in the square brackets does not have the Born pole O^^. The 
integration domain is shown in Fig. 2a. The upper bound on Qxk follows 
from the angular ordering. A further simplification results if we introduce 
logarithmic variables x = ln{P/k), z = \n{K/k), ( = In (t^/c/Qo), ^ = 
In (9ii-fcfc/Qo), A = \n{Qo/A), and the corresponding density p^^\x,C) = 
d'^n/dxdC = 2TTkH^p^p\k). Then, with a^{x) = (3'^/x 

pW{^^^C) = a\x + C + X)+ rdza\z + C + X) 1^ dip^^\z,i). (24) 

JO JO 
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This equation has a simple solution for constant as 0]. 

p«(fc)=A^p(fc)/o(2a/<). (25) 

The case of running as will be discussed in Section 3. 

This completes our first application of the pole approximation. The single 
particle density (pSD plays also an important role as the Green function or 
the resolvent for other equations, see Sect. 3.1. 



2.4.2 Simplified equation for connected correlation functions 

In the second example we use the pole approximation to reveal the simple 
structure and natural variables for the connected angular two body corre- 
lation functions. Momentum integrated connected correlations satisfy equa- 
tions analogous to Eqs. and ( pll) . Due to the different singularity struc- 
ture of the product term dprod{^i-,^2) and r*^^)(i7i, ^2) itself, the simplified 
integrals over the parent momentum K will be different in both equations. 



In (|2T|) the angular integral dVL^ is dominated by two regions K \\ ki and 
K II Choosing appropriate polar variables in both cases we get neglecting 
the nested contribution, 

^f{^i,^l2)= (26) 

^ f, fa\m,)p^^in,.) r [ej,,/^in,,j] + (i ^ 2). 

Vf JQq/^12 K Jkk ^Kki ^ ^ 

The singular part around each pole was integrated down to the elementary 
cutoff kk = Qo/K while the slowly varying function QKk2{^Kki,'&i2, is 
approximated by its central value i^u in all nonsingular expressions. Similarly 
QpK = 1,2. The upper bound for Qxki is chosen as ^12 since the 

contribution from beyond this scale is asymptotically negligible. 

It is readily seen from Eq.(PBD that, although d^"^^ contains a product 
of two poles (or more general, power singularities), the structure of A^^^ is 
simpler. The pole dominance in d^K integration splits A*^^-* into a sum of two 
terms where the parent is almost parallel to each of the final partons and 
consequently A*^^) is a sum of two contributions. It is also easy to see that 
A^^-* defined by the decomposition (|26| ) depends only on the two variables "di 
and ^912,"^ = 1,2. Hence they may be regarded as the natural variables for 
this process. The integration domains in this case are shown in Fig.2e,g. 
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This structure of the inhomogenous term A imphes simphfication of the 
parent integration in the integral equation for r^^-* itself , Eg . (pOD . It turns 
out that r^^) also decouples into a sum of two terms vf'^ (z = 1, 2) and each 
of them satisfies the following equation 

Again the pole dominance was used. The upper limit on Qku, follows from 
the angular ordering as in the single particle distribution Eq. p3|) , while the 
lower one is set by the minimal virtuality required for the subsequent de- 
cay into ki,k2 pair (see Fig.2f,h). This discussion provides a quantitative 
example of the general statement made at the beginning of Sect. 2. 3. Similar 
equations hold for higher order correlation functions (see Sect. 5). They all 
imply a remarkably simple singularity structure of the connected correlation 
functions. 

2.5 Generic equation 

At this point it is is worthwhile to emphasize a large degree of the universality 
in the description of seemingly different observables characterizing the QCD 
cascade. It turns out that many physically different processes are described 
by mathematically the same equation requiring only an appropriate renaming 
of variables and a proper choice of the inhomogenous term. For this reason 
we introduce a generic integral equation 

his, ^, P) = diS, ^,P)+ — —a\pT)h{5, VI/, K). (28) 

JQo/S K JS W 

where the scale of the running as depends on the kinematics of the specific 
problem. In principle one should distinguish two classes: a) pt = Kd{P'^), 
and b) pt = K'^. For example, class a) contains connected angular cor- 
relations, cumulant moments and doubly differential single parton density, 
and class b) correlations in the relative angle ^12 and single parton energy 
distribution. In fact equations describing both classes are related by simple 
integration, therefore it suffices to consider only one class. The inhomoge- 
nous term d depends usually only on two out of the three indicated variables. 
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It is essentially different for single parton spectra and for higher densities. 
Typically for angular correlations it has the following form for n-parton ob- 
servables in the running case 

d ~ exp (2n(3^\n{P6/A)^ . (29) 

Solutions of this equation are discussed in the Sect. 4. 



3 Multiplicity and single parton densities 

This chapter summarizes various techniques of calculating different single 
particle densities which we needed in the derivation of various multipar- 
ton correlations. Most of the results discussed here are known 0. The 
motivation here is to introduce our notation and the methods, which are 
used in Sect. 4 in the derivation of the angular observables. Furthermore we 
want to expose the emergence of the power behaviour in the variety of ob- 
servables. This aspect of the QCD cascade, i.e. the explicit proof of the 
self similarity /intermittency in the time-like partonic showers, has attracted 
considerable interest only recently [^|. Some results concerning the integral 



representations of various densities were not derived before. As we shall see, 
not all observables show selfsimilarity even for the constant . It is one 
of our goals to identify the conditions necessary for occurence of the power 
behaviour. 

Total multiplicity. Multiplicity of partons n{P, B) in a jet with momentum 
P and half opening angle is the simplest quantity which was studied in 
perturbative QCD |^. The equation which determines total multiplicity 

n(P,e)= I I ^p?(lnfc,ln^), (30) 

follows readily from Eq.(p3D upon integration over the parton phase space 



(|30D . Due to the logarithmic singularities of Born diagrams, it is convenient 
to use the logarithmic variables. Moreover, the structure of the inhomogenous 
Born term, together with the resulting equation, implies that total multiplic- 
ity depends only on the "virtuality" of a parent jet P9. One obtains, 

n{Xx) = h{Xx) + dv f dua'^{u + \)n{u + A) (31) 
JO Jo 
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X = ln(Pe/go),M = ln(K^/Qo),XA = X + A, A = In (go/A). (32) 

with the Born term b{Xx) = dv Jq dua^^u+X). For constant (a(x) — > a) 
one finds the simple solution 

n{X) = cosh (aX) - 1. (33) 

Hence in the high energy limit (P oo) multiplicity grows like a power of 
the virtuality 

" ^ (f )° 

with the anomalous dimension 70 = a. 

For running ag one solves the differential equation equivalent to Eq. (|32|) , 
which can be reduced to the Bessel equation by appropriate change of vari- 
ables. The solution satisfying initial conditions implicit in (p^) reads (H]) 



n{Xx) = 2(3^^Xx[Koiwo)h{2l3^JXx) + Ioiwo)Ki{2(3^ Xx)] - 1, (35) 

where wq = 2fi\f\. At high energy using /„(-2) ~ e^/ \J2t\z, Kn{z) — e'^^Jii /2z 
one recovers the well known expression 

'^<-'-^(^)"^^ (-) 

The results for constant can be found in the limit /3, A — 00 with 
/3 / = a kept fixed which will be referred as the " a-limit" . It also implies 

/ - v^e-^^^, a(Pe) - a, / (^^ j ^ (^^^ j (37) 

So in this limit Eqs. ( p5| , |36D reduce to (^,|^) respectively. It is seen that the 
power growth of the total multiplicity occurs in the constant as case. As 
is well known running violates this simple scaling and multiplicity grows 
slower than any power of the energy. 

Angular distributions. The angular distribution of partons in a (P, 6) jet 
is defined as 

p(^)(e,n)= r p^'\y,^,Y,)dy, (38) 

J — f 
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where Yx = Y + \, Y = In {P{}/Qo),y = In (k/Qo) and ^ = In {{}). Integrating 
Eq. (|23|) over parton momentum k gives the following equation for the angular 
density (m) 



= 6i(n + 0+ f dZa\Zx + i) f^dCp^'\CZx), (39) 



where Z = ln{K/QQ),( = ln(-?/') measure the momentum and the emission 
angle of the intermediate parent parton, bi{u) = P^ln {u/\). 

Eq.(^) implies that the angular density is independent of the jet opening 
angle G. This follows from the angular ordering, and the pole dominance. 
Namely all emissions from the intermediate parents (relative to the final par- 
ton k ) are limited by the final angle '& = Qpk between P and k. This applies 
to other distributions provided none of the angular variables is integrated. 
This simple property has rather important consequences, for example 

P«(e,n) = ^^^^^|.^,= P«(lA + 0, H = lne. (40) 

Therefore all results derived previously for the jet multiplicity can be directly 
translated for the angular density as well. In particular, the complete solution 
of Eq. (^) for running as reads 

Pp\^) = ^[i^o(^o)/o(2/?V^) - h{w^)Ko{2(3^x)]. (41) 

One may also derive this result solving the differential equation equivalent 
to Eq. (p9|) . In the high energy limit we get 




The corresponding exact and high energy results for constant are 



P?(^) = ^sinh(aln(^/«:)), P?(^) ^ ^ (^) (43) 



As for the total multiplicity, the power law emerges only at very high energies, 
and is violated by running as ■ Note that the angular density depends on Pd, 
i.e. on the transverse momentum of a whole jet with respect to the parton 
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direction. Result (^,^) is also relevant for the intermittency study. It shows 
that even the first factorial moment may have nontrivial intermittency index. 

Finally we quote the integral representation for running as which follows 
from the direct iteration of Eg. (p9D , cf. also Eq.(|77D . 

= + r dv r duR{Y -v,v -u,u + \)bi{u + \). (44) 



Where R{x, t, A) = dtR{x, t, A) can be calculated from the momentum distri- 
bution, cf. Eqs.(|8l|g|J^). 

Momentum distribution. Contrary to the angular density, the momentum 
distribution p^^\x,t,X) = dn/dx depends essentially on two variables, for 
which we choose x = In {P/k),t = In (Qk/Qo). As a result the corresponding 
integral equation, 

p(^) (x, t, A) = In + r dz f dr ^^—p^^^ (^, r. A) . (45) 



A JO Jo z + r + A 

has more complex structure, and its analytic solution is not known. It turns 
out, however, that the moments of the energy distribution can be calculated 
analytically. Consequently the distribution itself can be reconstructed. To 
this end define the dimensionless moments of the energy distribution 



Pn{X, A) = P^" / A; {k, P, e, X)dk = I e-'^V^^^ {x, X - x, X)dx, 
Jqo/& Jo 

(46) 

with X = In (PQ/Qo) and t = X — x. It follows from the master equation 
(^51) that pn satisfy the following differential equation 

/i„(X) + n-^/i„(X) = -^(p^(X) + 1), (47) 



dX^"^"' ' dX^"' ' X + A 
with the boundary conditions 

/in(0)=0, A^„(0) = 0, (48) 

Eq.(P7|) is easily transformed into the confluent hypergeometric equation. 
The solution satisfying the boundary conditions (|48|) reads. 



/i„(X,A) + 1 = 

Y[a)Ze-^ [$(a - 1, 1, nA)F(a, 2, Z) + F(a - 1, 1, n\)^[a, 2, Z)\ (49) 
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where Z = nXx, a = l + /?^/n, and F and $ denote the regular and irregular 
confluent hypergeometric functions respectively. Inverse Laplace transform 
of® 

p^^\x,X -x,X)[e{X -x) -e{x)]= /^e"Xs,X,A), (50) 

can now be used to derive the asymptotic behaviour of the momentum 
spectrum at high energies. To this end we perform the s integration in 
Eq. (|50|) by the saddle point method. The saddle point s* {X, x) — >• for 
X — > oo,x/X = const., therefore it is sufficient to use the appropriate 
asymptotic forms of the confluent hypergeometric functions (see Appendix 
A). Using Eqs.(|234],|235D we flnd the following condition 



X , S _ . X , r, ,9n/1 1 

1 n , 4ag , 




(51) 

A = s"^ + (3^/{X + X) , 5 = + 4ao - 2s/A (52) 

with Oq = /3^/A. This determines the saddle point s = s*{X,x). The integral 
( |50|) is dominated by the neighbourhood of this point with the following 
leading result 

p«(a;,X-x,A) ~ exp |^2s*a; + 

exp(^(v/A-v^)(l-^:^!±|I)) exp 

This is the double logarithmic expression for the famous hump-back plateau 
which was conflrmed experimentally and is considered as an important test 
of the perturbative QCD [|15] . 




We conclude this subsection by pointing out some interesting limiting 
cases of the above expressions. 

1. Total multiplicity. The zeroth moment of (^) gives the total multiplicity 
of a jet. Indeed using the asymptotic form of the confluent functions p6| 



lima-,oo,b=constF{a,b, x/a) = r(6)x(^ ^'^^'^Ib^l{2^/x) and lima-.oo,b=const[^{l + 
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a — b)^{a,b,x/a)] = 2x^^~''^^'^Kb^i{2y/x), one readily obtains the already 
quoted result (|35|) for n. 

2. Constant as limit. This limit, the a— limit cf. Eg. (P7D , is realized formally 
with A,/9 — > oo,/9^/A = = const. After some algebra one obtains from 
( ^91) much simpler expression (see Appendix B). 

/x„(X) + l = 

1 ^ ^ ^ ^ ^ X 



(V + 4a^ + n) exp [ — {y/n? + 4a? — n)] 



+ (Vn2 + 4a2-n)exp[-y(Vn2 + 4a2 + n)] ) (54) 

which can be also derived from the Laplace transform of the analytic solution 
( 1^ or by solving Eq.(^7D for the constant . 

Finally we consider the a— limit of the hump-back formula (p4|). Equation 
( |5TD simplifies considerably and one readily finds 

= - . (55) 

Jx{X — x) 

and the asymptotic form of the density (p3) reads 



p(^)(a;,X-a:) ~exp(2a^x(X-x)). (56) 

This agrees with the asymptotic form of the exact solution for the constant 
Us which follows for example from Eq.( 



p{x,X-x) = - l/i(2ai/x(X-x)). (57) 

We conclude this Section with two observations. 

1. The momentum distribution, Eq.(|56D does not show the power behaviour 
in any of the two variables. Therefore not all observables (and processes) 
are revealing the fractal nature of the QCD cascade. The selfsimilarity is 
destroyed by the second (in addition to P) independent scale k which exists 
in this process. This observation suggests, that only processes containing a 
single scale are suitable for searching for the intermittency/selfsimilarity. 

2. Note also that if the second scale is fixed to be proportional to the first one 
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X = eX say, than we recover the power behaviour. In our example Eq.(|5 
becomes 

p«(x,X-x)~ (^^^j , (58) 

and the intermittency indices depend on the ratio e = x/X. We shall find 
this phenomenon also in other more complicated situations with many scales. 



4 Methods of the solution for higher densities 

4.1 Solutions in terms of the resolvents 
4.1.1 General resolvent representation 



There exists a rather general recursive scheme for solving Eq. (|I5|) for a 
density of arbitrary order. First, observe that the inhomogenous term dp'\ 
which corresponds to the direct emission from the parent parton, is built 
from various products of the correlation functions of lower order. Only if all 
n derivatives act on the Z function under the d^K integral in the exponent 
of Eq.(|^ can we recover the n-th order density. This gives the last term in 
Eq.(|l5|). In all other cases the n derivatives are split among various factors 
resulting from the differentiation of the exponent giving various products of 
the densities of lower order. Integral of any particular density, of lower than 
n-th order, can be always replaced by the density itself and corresponding 
inhomogenous term of yet lower order, by using appropriate integral equation 
for the lower order. 

Second, the unknown function p^"'^ appears only in the last term in 
Eq . ([15|) . The simplicity of this structure allows us immediately to solve 
equation (|1^) by iterations 

Pp^iJ^li ■■1 ^n) = 

CO „ 

/ d^K,.J^KrMp{K^)MKAK2)--MK^_,{Kr)d'^;^l{h,..,K), (59) 
or symbolically 

p(") = (-^)orfW. (60) 
1- M 
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Further, we note that all but one Km integrations (m = l,...,r — 1) can actu- 
ally be done, since the complete Km dependences are given by the Born cross 
sections and variable boundaries of the inner integrations. Hence changing 
the orders of integrations we get the following integral representation of the 
solution 

pPih, .., kn) = dPih, .., kn) + Rp{K, a)dP{k,, .., kn)d^K, (61) 
where the resolvent Rp{K, a) is given by 

oo „ 

Rp{K,a) = Y^ / d'Kr.i...d'K^Mp{K^)...MK^_,{K). (62) 

r=l "'r,0x^_jx>o- 

For constant ag, Rp{K,<j) can be calculated exphcitly (see Sect. 3 for more 
details) but the results ( |6T| , |62D are valid for running as well. In these 
formulae a is the minimal opening angle of a jet K and has to be de- 
termined for each case separately; in general it depends on all momenta 
a = a{ki, . . . ,kn, K), for example, a = QkiK for n = 1. 

The role of a can be better understood if we consider in detail the con- 
tribution from the second iteration, see also Fig.4, 

pp'\}C)=f d^'K.MpiK,) I d'K2MKAK2)d^Kli}^). (63) 

Fi and F2 denote appropriate boundaries resulting from all three constraints 
Vp{Ki), Vki{K2) and F7^2(/C), where /C denotes collectively all momenta of 
the final partons /C = {/ci, . . . , kn\- After changing the orders of the Ki and 
K2 integration we get 



p(,"'2)(/c) = / SK2 L d'K,Mp{K,)MKAK2) 



41 (/C), (64) 



which allows to identify the term in the square brackets with the second 
order contribution to the resolvent, Eq. (p2|). However, because the Ki 



integration is now done at fixed K2, its boundaries are influenced by the 
constraints implied on K2 by the configuration of the final partons /C. In 
particular, angular ordering requires that QkiK2 > c"(^) -^2) with cr defined 
as above. Similar considerations lead to the general formula, Eqs. (|6^j68| ), 
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with a being the lower cutoff of the emission angle QpK and, at the same 
time, the opening angle or a measure of the "virtuality" of the jet K. 

Eq. ( |BTD has a simple interpretation in terms of the cascading process 
(see Fig.l). In fact, the resolvent Rp{K,a) is nothing but the inclusive 
distribution of a jet {K,a) in a jet (-P, 0). Note the dual interpretation of 
the angles a and 9, they denote opening angle of the corresponding jets, but 
also they are simply related to the viruality of the parent partons. 

Eqs.(^j61|) give the recursive prescription for explicit calculation of general 
(fully or partly differential) multiparton densities of arbitrary order. 



4.1.2 Resolvents for fixed as 

The explicit form of the resolvent can be readily derived by considering the 
integral equation for R itself. It can be written in many equivalent forms. 
We shall use the forward master equation for the density in the logarithmic 
variables a; = In (P/ K),t = In (■(^/cr) . 

P(x,t) = ^ = ^^^ (65) 
dxdt Mp{K) 

From Eq. (|62D we have (symbolically) 

R = M+MoR, (66) 

or explicitly 

R{x,t) = + r dz f dTR{z,T)., (67) 
Jo Jo 

The pole approximation was used similarly to Sect. 2.4.1 to simplify the 
three-dimensional integrals. Ka is the virtuality of the {K, a) jet emitted at 
the polar angle ^. Equation ( p7D is identical with Eq.(p4D for constant 
with the replacement Qq Ka. Hence the solution (|68|) follows. 



Rp{K, a) = Mp{K)Io{2a^HP/K) ln(eKp/cr)), (68) 

Solution ( |6T[ ) together with the formula (p8| ) provides useful integral ex- 
pressions of many observables for constant . We quote here several re- 
sults leaving other more involved applications for subsequent sections {Y = 
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ln(P/Qo),l/ = ln(A;/Qo),S = In 6, ^ = In ^9). 

1. Total jet multiplicity and total factorial moments. 

n{X) = b{X) + a sinh {a{X - u))b{u), b{x) = -a^a;^ (69) 
Jo 2 

f„{Y + S) = + S) + £^ dCIoi2a^{Y-Z){E-C))bn{Z + (), 

(70) 

where fi{X) = n{X),bi{X) = 0^X72,62 = exp(2aX), and higher are ob- 
tained recursively, (see Appendix D). 

2. Energy density p^^\y, Y, 5) = dn/dy. 

p^'\y,Y,E) = biE + y)+ ^ dZ T ciao(V(^ - - C))% + 0, (71) 
with 6(X) = a^X. 

3. Angular distribution p^^\^,Y) — dn/d^. 



pW{^,Y)^b{Y + + dZ dQh{^{Y -Z){i- Q)b{Z + (), (72) 



where h{X) = a^X. 

4. Inclusive density in both variables p^^\y,^,Y) = dn/dyd^. 



P^'\y, Y)^a'+ r dZ f dCIo{^{Y-Z){^-C))a'. (73) 
Jy J-y 

5. Energy-energy correlations p^'^\yi,y2,Y,'E) = d^n/dyidy2- 



p^'\yr,y2,Y,E) ^b{yi,y2,Y,E)+ [ dZ T dCIo{2a^{Y -Z){E- (My,, 2/2, Z, Q, 

Jy> J-y< 

(74) 

where 

b = pW(yi,r,s)pW(i/2,r,s) + '''I^^^y^p 2{2a^\yi-y2m + y<)), (75) 



pW(y,y,S) = a J ||±^/i(2aV(y-y)(S + y)) 



(76) 
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and y>{<) denotes the greater (smaller) one of the two "rapidities" yi,y2- 
Some of these expressions can be further simplified, for example the explicit 
forms of the total multiplicity and of the single parton densities quoted in 
points ( 2 - 4 ) are well known 0|. 

Equations (|7l] - ^ ) can be derived, either as the resolvent solutions of 
the appropriate integral equations, or as the pole approximation of Eq. (|6T|) . 
For example, solution (|6l|) of Eq. (^) for the single parton density reads 

p?{k) = Mp{k) + J d'KRp{K, eKk)MK{k). (77) 

In this case the minimal opening angle of a {K, a) jet a = Qku- The integral 
is saturated by the 0^-^ ~ pole of the inner bremsstrahlung. Approximat- 
ing QpK\Qj^^^Q ~ 9pfe at the pole, and changing variables gives Eq.([73D. 

4.1.3 Resolvents for running as 

The forward master equation for the resolvent reads in this case 

There are two novel features: the coupling constant a? depends on the ap- 
propriate transverse momentum, and the resolvent depends on the three di- 
mensionless ratios as opposed to the constant as case. In the logarithmic 
variables x' = In {P/K),t' = In {'d/a), X' = In (Ka/A) Eq.(^) takes the form 

R{x', t', A') = a\t' + A') + r dz t dra^z + t' + X')R{z, r. A'). (79) 

Jo Jq 

It will prove useful to introduce the integral over the polar angle R{x', t', A') = 
/q dTR{x' , T, A') which satisfies the equation 

R{x', t', A') = In + r' dz f dr' J' .M ^, r'. A'). (80) 

A Jo Jo Z -\- T -\- A 

This is equivalent with the integral equation for the inclusive momentum 
distribution of partons (|45|). Therefore we conclude that 

R{x',t',X') = p{x',t',X'). (81) 



22 



It is worth to emphasize the different physical meaning of the parameters A 
and A'. While in p: A = ln(Qo/A), on the LHS: A' = ln(i^a/A). This is 
a natural and important generalization from the elementary partons to the 
virtual jets. While for partons Qq is the fixed cut-off, for jets its role is taken 
by the lower scale (virtuality) from which we are following their perturbative 
evolution. 

To summarize, we find the resolvent R{x,t, A) by solving Eq.(^) for the 
inclusive momentum distribution of partons with the appropriate reinterpre- 
tation of variables. This will be further discussed in the next section. 

Once the resolvent R{x',t',X') is known, any observable can be directly 
calculated from the representations analogous to Eqs. (ffll - [73| ) also in the 
running ag case. 

4.1.4 Saddle point solution of the generic equation 

We are now ready to solve the generic equation (^). Only class b) be will 
be considered. First, observe that the equation 

rP dK n 

hn{5,i},P) = d45,P)+ — —a'{K^)h^{5,^,K), (82) 

jQo/S K Js W 

is identical with the the integral equation for the energy distribution 

p{k, P, 6) = b{k, 9) + / — / —a'impik, K, M/), (83) 

Jk K JQo/k W 



if we replace in (83) 



k ^ 5,P ^ ^ P,^ ^ K,K -^^,b^ dn. (84) 
Together with Eq.([77D this gives the following resolvent solution of Eq.(^2D 

fP dK r^d^—f'd P K'^\ d 

hjs, n ^ us. p) + - / -fl -. -) _d„(^. Ky 

(85) 

The leading contribution is yet simpler (see Appendix C for the details) 
, ON , (P^\ dK-(d P K5\ ^ (K5\ 
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Using (^) for the energy moments we get from the inverse Laplace transform, 
Eq.(|50|) and from Eq.(0) 

with A{s, A) = r(a)$(a - 1, 1, sA), /(s, X) = sXe-'^F{a, 2, sX) and a as 
in Eq. (|i9|). 0. Performing both integrals in the saddle point approximation 
gives our final result (see Appendix A) 



/i„(5, P) = VnVn exp [2P^\n {P^/A)uj{e, n)j , (88) 

where 



_ In {^/6) 
^ ~ ln(P^9/A)' 

The universal function Ci;(e, n) is given by 



(89) 



ijj{e,n) = 'y{z{e,n)) + ez{e,n), (90) 



1 



7(^) = -(vV+I-z), (91) 

and z{e, n) is the solution of the following algebraic equation 

'^^{z)-2\n-i{z)-ez'^ = n^ -2\nn. (92) 

Equations (pO|-|9^) completely determine the scaling function c<j, in particular 
its behaviour around the end points of the kinematic range can be readily 
obtained 

w(e,n)~ ri-^e, 6-0, (93) 



;e,n)~ v^r^y^ln + n2 - 2 - 2 Inn, e ~ 1. (94) 



In Fig. 5 we show the results for uj{e,n)/n with n = 2,3,4. They approach 
a finite limit for large n, e fixed, see Eq. (|110|) below. The normalization of 



^The contribution from the second term in Eq.(^9|) is negligible for large X. 
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( pHD can be conveniently split into two factors. The first one is independent 
of the inhomogenous term 



V 



n 




l\ 7(z)2 - I + ez{2-f{z) + z) 



(95) 



while the second one Vn represents the normalization of the inhomogenous 
term dn , hence it depends on the considered process. It can be readily 
calculated for the cases of interest (cf. Chapter 6 and Appendix A). At fixed 



e both factors provide a nonleading (of the relative order 1/yln (P-i^/A)) 
corrections to the asymptotic scaling function uj{e,n). Since Vn{0) = 1 and 
it is regular at e = it remains close to one for small e. V2 rises up to 1.6 
for e = 0.5 and Vn < Vn-i- 

4.2 WKB solution of the corresponding 
differential equation 

4.2.1 The case of running as 

An alternative method to solve the generic equation (p8[ ) consists in solving 
the corresponding differential equation (for pt = K"^) 



in the logarithmic variables u = ln{P6/A) and v = ln{'i}/6) with boundary 
conditions along u = ln((5o/A) and v = 0. Whereas we obtained from the 
resolvent method the exact asymptotic solution, the WKB solution yields 
easily power expansions and nonasymptotic corrections. 

As we are interested in the exponential (power) behaviour and corrections 
to it, we write a differential equation for the exponent and study its expansion 
in appropriate variables (generalized WKB-method) 




S{u, v) = In h{u, v) 



(97) 



5^S 




(98) 



SuSv 
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with = fP'jiu + v). The boundary conditions are obtained from the 
inhomogenous term (pO]) as 

^(m, 0) = 2n/3v^ (99) 

In our high energy approximation the second boundary conditons for u = 
ln((5o/A) will be satisfied only approximately. 

In order to find convenient expansion variables for the solution of (03) 



with ( p9D we perform the first few iteration steps of the generic equation (^) 
with (129|) . In the transformed variables X = 2ni3y/u = 2n/3 ln^''^(P(5/A) and 
V = ln(79/(5) and after expansion of a^(if\E') this can be written in the form 




m 



h{X, t^) = + 2/32 / — / ^ c„ ( — j h{x, v')dv' (100) 

Starting with h^i^X^ v) = exp X one finds for the first iteration after successive 
partial integrations in X 

2 

hi{X,v) = e^ — (coo + cio^ + C2o^ + ...) 

2 

+ e^^(coi + Cn^ + C2i^ + ...) + ••• (101) 

Here we neglected the contribution from the lower limit Xq of the in- 
tegral, i.e. we require u (or the momentum P) sufficiently large. In every 
higher iteration step for /i„ the factor v/X = X ■ [v/X"^) multiplies hn-i- 
The expressions for /i„ like ( |101| ) suggest a double expansion of h in l/|3^/u 



and v/u. Instead of calculating the coefficients directly from the iteration 
we expand the eikonal S in these variables and determine the coefficients 



from the differential equation (|98|). However v/u is not a good expansion 
variable as it can vary between and 00. Therefore it is more convenient, 
but equivalent, to replace u ^ u + v and to choose as expansion parameters 

" - : „/ - . ' (102) 



u + v ln(P?9/A) (3Vu + v [3^\n{P^/A\ 

with < e < 1. This choice also leads to a full cancellation of the a?{u,v) 
term in Eq.(p8|) in leading order of C/^ ^ {u + v)~^ and e and to a rapid 
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convergence of the expansion in the region e < 1/2. Therefore we write 



1 oo oo 

^(^'0 = tEEc^A^' (103) 

j=ofc=0 

r\ OO 

= 7E^.(^)C^' (104) 

The coefficients Ckj can be obtained successively from (p8|) ; the boundary 
condition (^) imphes a;o(0) = n. Whereas the expansion of h in ( |101D would 
involve also negative powers of C, the limitation j > in ( |104|) is required by 
the boundary condition for small u. As u itself should be sufficiently large 
we consider the limit ti/w-Clore~l. Consistency with ( ^9]) can only be 
obtained for nonnegative j. In the high energy limit C with e fixed we 
obtain 

S{u,v)^ 2u^{e)/C (105) 

Note that this limiting form agrees with the result from the resolvent method 
Is]). Inserting the expansion (|104|) into (|98|) we can also obtain differential 



equations for ujj{e). 

In leading order O(C^) the double differential first term in ( p8l) drops and 
one finds for uj[e) = ujois) the WKB equation 

{uj-2euj'){iu + 2{l-e)uj') = 1 (106) 

which easily yields the Maclaurin-expansion for u 

-(^.") = « - (1 + + + + 3)^^ + ■••). (107) 

The differential equation, Eq.( 106 ) together with the boundary condition, 
corresponds to our algebraic result Eq. ( ]90| - |92D . As a consistency check we 
have reproduced the first four terms of the expansion ( |107|) starting from the 
definition (^01-0). 

For the terms nonleading in ( we obtain in lowest e orders 

/ , 1 - 3 2 

1 3 
^2ie,n) = + u^ie^n) = —e + ■ ■ ■ . (108) 
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Alternatively we may solve Eq.( p.06| ) in an expansion in 1/n; note that 
Eq. (|107|) suggests u;{e,n) = nf{e) +g{e)/n. The leading term n/(e) is found 
as nv^l — 6 and solves Eq.( [106| ) if the r.h.s. is neglected. Then one obtains 
for g{e) 

g'{e) + g{e)/{2 - 2e) - l/{2Vl^e) = (109) 
with (7(0) = 0. Finally one finds 



= nva^ (l-—ln(l -£)) + •••. (110) 

Remarkably, this approximation has an accuracy of better than 1% for e < 
0.95 already for n = 2 and is therefore well suited for numerical analysis. 

The same type of asymptotic behaviour ( |105|) was found for the cor- 
relations studied in |]2T[], their effective power u^ff is related to our u by 



^eff = {n- w(e, n))/ (1 - VI - e)- 

So far we have analysed the generic equation (^) with pt = K"^ which 
applies exactly in this way to the 2-particle correlation p^'^\'&i2)- For the 
cumulant correlation T^"'\{^li}) the same equation holds but with pt = K'd; 
also the inhomogeneous term obtains a variable prefactor 

A{u, w) ~ (m + vy^u'^exp{2n(3^) (111) 

Taking into account the mixed scale Kd and writing F = exp 5, A = exp S' 
we arrive at the modified differential equation 

dSdS dS 1 dS 
_ ^ \ \ 

du dv u + V dudv u + v du 

^s-s' f 9S' ^ dS'dS' ^ 1 dA\ 
\ dudv du dv u + v du j 

The form (|TTT]) for A and ( TOITTUD for S implies S - S' h'^{e)^/^rTT^ 



for large u + f , therefore an exponential decrease of the r.h.s. of ( |112| ). In 
leading order of = [u + v)~^ only the first two terms of the left-hand- 
side contribute as before, therefore one obtains the same asymptotic solution 
( |105| ), and differences appear only at the nonleading level. 
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4.2.2 Results for fixed a 



We start from the generic equation (pS]) with constant a but with inhomoge- 
neous term 

c/~exp(naln(P(5/go)) (113) 

The approriate expansion variables are now u = ln(P(5/Qo) and v = ln{'d/6) 
as before, also e = v/{u + v). As above we arrive at the same differential 
equation (^) but now with boundary condition 

S{u, 0) = nau. (114) 

Similar reasoning as above leads to the ansatz S = a{u + v)uj{e,n) where 
uj{e,n) satisfies a differential equation as in ( |106| ) but without the factors 2. 
The solution is found as Lj{e, n) = n — {n — ^)e, so finally 

S = an{u + v) {l - {1 - n-^)e) (115) 
= a {nu + n~^v^ (116) 

As is easily seen this linear function in e is actually the exact solution of 



Eq.(p8|) with ( |114| ), i.e. of the generic integral equation ([28| ) where the con- 
tribution of the lower limit from the K integral is neglected. So the running 
of the coupling is most directly seen in the nonlinear behaviour of the uj{e, n) 
function. 



5 From integrated to differential correlations 
5.1 Multiplicity correlator Fg 

The results in this subsection are needed later for consistency checks of our 
differential results; they are therefore rederived using our methods although 
they are largely known already (see @], for example). 

The second factorial moment, or the total number of pairs in a jet, pro- 
vides the normalization of the two parton densities 

/2 =< n{n - 1) >= J d^kid^k2P^p\ki, h) (117) 
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Its integral equation follows from Eq. (JT^) for the two parton differential den- 
sity, X = In (Pe/go), = X + A = In (PB/A) 



f2{Xx) = g2{Xx) + r dv r dua\u + \)h{u + A). (118) 
JO Jo 

This differs from Eq.(^) for the total multiplicity only by the form of the 
inhomogenous term g2 = J d^kid^k2dp\ki, with d = dprod + dnest given 
by (^[12]). For constant as one gets exphcitly g2iX) = rf + 2n — aX^ and 
the solution of Eq. (|118|) reads 

/2(X) = ^(cosh(aX) - 1)^ ^ 1 (^^^ . (119) 

The normalized correlator F2 = /2/'^^ is then found from ( P^ ) as P2 = | 

Using the resolvent representation Eq.([70|) one can easily show that for 
constant 

/,(X) = P, < n(X) X^oo, (120) 

with < n >= exp(aX)/2 and the coefficients Fg are given by the recursive 
relation (see Appendix D for the details) 



^.=7^Er. ^>i^ ^1=1- (121) 




Equation ( |120|) generalizes to the case of running as ■ In fact, the inde- 
pendence of the coefficients Fq in ( |1 20| ) of the running coupling (or more 
precisely of A and P) is sufficient to assert that the recursion (|121|) , derived 
for constant , must also hold for the running . Otherwise the usual 
correspondence between the constant and running coupling results would be 
violated. 



5.2 Correlations in the relative angle 
5.2.1 Definitions and integral equation 

The distribution of the relative angle between two partons in a jet with 
momentum P inside a cone of half opening angle B is defined as follows (see 
Fig.6) 

p(2)(^,2,P,e) = J pf{k^,k2)5{Qk,k,-^i2)d^kA. (122) 
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This quantity shares the simphcity (small number of variables) of global 
observables and, at the same time, provides the first nontrivial information 
about the differential structure of the correlations. We also consider two 
normalized densities. Firstly, 



4?„,(^i2,p,e) 

where the normalizing quantity in the denominator is defined as in ( |122| ) but 
with Pp\ki,k2) under the integral replaced by Pp\ki)pp{k2); then, for an 
uncorrelated distribution p^'^\ the normalised correlation is r{'di2) = 10. 
Secondly we define 

with the full multiplicity in the forward cone as normalizing factor. The 
importance of these different normalizations will become clear in Sect. 7 
where f is found to exhibit better scaling properties. 

The correlation function is defined in (|122|) as an integral over the full 
forward cone, therefore it matters in general whether the leading particle is 
included or not. We will derive first the correlation without leading parti- 
cle p{^i2), the one with leading particle pi'&u) is then obtained from Eqs. 
(|l|,|mD as 

P^p\A2) = pf{^i2) + 2p2)(^i2). (125) 

Whereas these two cases are clearly separated in the theoretical analysis it 
is not so straightforward to realise experimentally the case with the leading 
particle subtracted. One could think of a flavour tag, for example subtracting 
the charmed particle in a charm quark jet. 

For clarity we consider first the gluon jet and come back to the quark jet 
at the end of this subsection. 

The integral equation for p^p\^i2) is obtained by integrating Eq.(|T^ over 

■^This type of normalisation is also suggested in ("correlation integrals"). Exper- 
imentally this normalisation is obtained conveniently by "event mixing" as is often done 
for the determination of Bose-Einstein correlations. 
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the phase space of final partons at fixed i!}i2 

(126) 

The derivation of this equation is similar to that of Eq.(^). Now however the 
polar angle = Qpk and as a consequence the relative angle {^u determines 
the lower bounds for \l/ and K integrations in accordance with general rules 
discussed in Sect. 2.4.1. Namely, the leading logarithms always emerge 
from K II ^1(^2)- In this configuration the minimal virtuality of a jet K, 
which subsequently emits a pair {ki,k2) with given angular separation, is 
controlled by the relative angle -t^u- Note that the individual directions 
fifcj, are already integrated over, therefore the "most narrow" singularities 
©i^Tfci, 0i<:fc2 ~ Qo/K are included in p'^'^\'di2,'^ ,K). One can say that the 
relative angle di2 is the only remaining scale controling the singularities of the 
integrand in Eq.( |126| ). On the contrary, for more differential distributions the 
boundaries depend in general on other different scales entering the problem. 
The direct term in Eq.( p.26| ) is defined in analogy to ( |122| ) 

d(^)(^i2,i^,e) = j df{h,k2)5{Qk,k,-^i2)d'h(fk2. (127) 

where dp {ki, ^2) is given explicitly by the inhomogenous part of Eq.(|l^) and 
can be written as d^^^ = d'^pj^^ + d^lst as discussed in Sect. 2.2. They can be 
derived from the angular distribution pp\^) 

d'^L{^u,P,e) ^ p?iA2) r'dAp?{A) + il--2), (128) 

J K 

rfil(^i2,P,e) = r fip£^(t9i2) r ^a^(Mi), (129) 

where k = Qo/P. The integration of the product term over the angles has 
been done in the pole approximation as explained in Sect. 2.4.2, where the 
first term corresponds to the pole for = 0, ^2 ^ "^u- 
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5.2.2 Exact results for constant as 

For constant the solution of Eq. ( |126| ) can be written as 

(130) 

with the resolvent 

J(x, y) = I,i2a^ln{x)ln{y)). (131) 

If the leading particle is not included the direct term is easily calculated from 
Eqs. (|128| , |129|) with p^p\'d) from (H) 

4J„,(79i2,P,e) ^ |^sinh(aln^){cosh(aln^)-l}. (132) 
c^il(^i2,P,e) = ^{cosh(aln^)-l}ln^. (133) 

(2) 

Note that integrating the direct term dp^od o'^^r the whole range of the relative 
angles gives the square of the total multiplicity in a cone 

r dA2d'pJUA2,P,&) = (cosh(aln(Pe/Qo)) - 1)' = n(Pe)^ (134) 

as it should. This shows that the pole approximation correctly identifies all 
leading logarithms also in this case. 

With the direct term given by Eqs.( |132| ) and ( |133| ) the integral in Eq.( p.30| ) 
can be done analytically (see Appendix E) and our final result without leading 
particle reads 

p^'\A2,P,e) = ^ E l/^'-^'WW - ^sinh(aln^), (135) 

^12 ^0 ^12 1^ 

where y = 2^\n (^^is/zt)/ In (6^2), z = 2a^\n {^2/^) In (6^2), and Ii{z) 
is the modified Bessel function of the order /. The same result (\Vd5\) can 



be obtained also by solving the partial differential equation corresponding 
to Eq.( |126| ) which is separable in the variables y and z together with the 
appropriate boundary conditions. 

If the leading particle is included one has to drop the (-1) in Eqs. ( |132D - 
( |134| ) and also the second term in Eq. ( P-35| ) as is found from relation ( |125| ). 
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It is easy to check the backward compatibihty of these resuhs with the 
second factorial moment /2 discussed in the previous Section. Integrating 



Eq.( |126D over the relative angle one readily recovers the integral equation 
(|118|) for /2. Also the direct integration of our solution ( |135|) gives back the 
expression (|119|) for the total number of pairs (see Appendix F ). 



In the high energy limit, k = Qq/P — >■ with kept fixed, the r.h.s 
of Eq.( [L35|) exhibits a power behaviour irrespective of the leading particle 
contribution, (see Appendix G for the details) 

which proves a selfsimilarity property of the QCD cascade for fixed as- Note 
that this simple behaviour emerges only for the well developed cascade, i.e. 
for sufficiently large angles 'di2 3> k. Eq.( |136| ) represents the two components 
of Eq. (|130| ): the first two factors correspond to the direct term {d ~ dprodji 
the last factor to the enhancement from the emissions of the intermediate 



parent jets. The second term in our solution (135) which distinguishes be- 
tween including or not including the leading particle is nonleading at high 
energies ~ {'&i2li^Y- The asymptotic result ( |136| ) can be obtained also eas- 



ily by solving the integral equation (|126|) , for example by iteration, with 

y(2) 
prod 



the asymptotic form for the direct term d*^^) ~ '^irid — {(^/'^'^i2){'&i2/ t^Y"', 



neglecting the nonleading nested term ~ {^12/ i^Y- 

For the normalised density r{^i2) in Eq. (|123|) the leading power from d^^)^^ 
cancels and one obtains 

/ e \t 

r(^i2)^ h^) (137) 



which is sensitive only to the nonleading exponent a/2. The leading exponent 



is not canceled for the alternative normalization Eq. ( P-24| ) in which case we 
find 



3a 



f(^i2)~— — . (138) 



^12 V^12/ 

It is noteworthy that these results are infrared safe as they don't depend on 
the cutoff Qq. 
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5.2.3 High energy behaviour for running as 

For running we derive the asymptotic behaviour of the density ( |122| ) using 
methods discussed in Sect. 4. The integral equation (|126|) is indeed one of the 
special cases of the generic equation with the substitutions 

6 ^^12, ^^e, h-^'di2pf (139) 

The inhomogeneous term is calculated from the one particle angular distri- 
bution in high energy approximation Eq.(^2|) using pole dominance Eq.( p.28| ) 
and neglecting the nonleading nested term as before. We obtain 

^ ^ exp (^4/?yMAWA)) (140) 



which has the form ( pOj) considered for the generic equation. Hence, applying 
the results of Sect. 4, the high energy asymptotic behaviour is governed by 
the function a;(e, 2) with the scaling variable 

e = (141) 
ln(Pe/A) ^ ' 

and we obtain for the unnormalised and normalised correlation functions 

P^'\^i2) = ;r^exp (2(3uj{e,2)J\n{Pe/A)) (142) 



2^9i2 V / 
r{A2) = exp (^2[3^ln{PQ/A) (u{e, 2) - 2VT^)^ (143) 

In the normalized correlation r{{}i2) the leading term 2\/l — e in the exponent 
cancels in the approximation (|11CI|) and we find in this case 

f3 



r{^i2) ^ exp |^-|^ln(P0/A)Vr^ln(l - e)j . (144) 

As to the correlation f defined in Eq. ( |124| ) we note that it behaves like 
f'i^u) ~ ^/"^u, therefore it is more convenient to consider the distribution 
differential in the logarithmic variable e, i.e. f(e) = 't9i2P*^^'*('(9i2) ln(PO/A). 
Using Eg . (|36D for n we find 
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f(e) = 2l3^\n{PQ/K) exp (^2/3^1n(P0/A) (cj(e, 2) - 2) j . (145) 

The prefactors in Eqs. (|142| , |143| ) and (|145|) provide the proper normahzation 
at e = consistent with ( |14CI| ). We stress that they are nonleading and 



therefore less rehable. Nonleading terms beyond DLA are not included in 
the present approach. 

For sufficiently large angles t9i2 (small e) we apply the linear approxima- 
tion a; ~ 2 — 3e/4 and obtain 



f2 /pp)x 4a{Pe) . . -|a{P0) 

for the normalized correlations r and f we recover the power behaviour in 



the relative angle as in the case of fixed as, Eqs. (|137| , |138|) , but with expo- 



nents which depend on energy and angle through the anomalous dimension 



a{PQ) = \jQas{PQ)/'K- The formulae ( |136| - |138D for constant are recov- 
ered in the a-limit A,/3 ^ cxd (e — *• 0), see Eg . (|37|) . 

It is worth pointing out the different scaling behaviour for the cases of 
fixed and running a^. In the latter case a scaling limit is approached in e for 
the quantity 

/ln(Pe/A) ^ 



2/3(cj(e,2) -2v/r^) . (147) 



For fixed cts we define 

_ In (6^2) 



^ in(pe/go) ^^^^^ 

i.e. replace A by Qq in ( |141| ); then in the limit e fixed, P — oo we find again 
the result ( p.36| ) or 

ln(Pe/Qo) 2^'- ^'^'^ 
The comparison of both cases shows (neglecting the difference between e and 
e) that ffist, the scaling behaviour occurs for differently normalized quanti- 
ties, and second, the limiting functions are different: the linear behaviour 
in (|147|) (namely \f3e) is approached only for small e. Therefore the study 
of the quantity r{di2) allows for some novel tests of the running of as in 
multiparticle final states and the same is true for the correlation f{{}i2) in 
complete analogy. 
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5.2.4 The behaviour for small relative angles 

If the relative angle approaches the lower cutoff n = Qo/P the DLA 
becomes more and more unreliable. Here we consider the limit of small 
angles ^^12 — ^ 1^12 ^ © but with ^912 still sufficiently above the cutoff k. 
In this limit the leading particle plays an important role, in particular for 
the correlation r{i)i2) which is considered ffist. 

For constant ag this limit in Eqs. ( p.35|J132|) corresponds to the lowest order 



in the coupling a , i.e. to the Born term. Without leading particle one obtains 

P?\^i2) ^ f/il(^i2)^(a7^9i2)ln2(^i2/«:)ln(eM2) (150) 
4™.(^i2) ^ (aV^i2)ln^(^i2/«:). (151) 

Both distributions vanish for ■&12 — > k but the nested term dominates by one 
power in In(?9i2/K). Then the normalised correlation becomes singular in this 
limit r(?9i2) \n{e/ k) / ln{'&i2/ n) or 

K\n(Q Ik) 

r{i!^i2) — ^ (without leading particle). (152) 

V12 - 

On the other hand, if the leading particle is included the additional 
term 2p^p\'di2) ^ (2a^ /^i2)\^{^i2/ ^^) in ([1251) dominates both pf^ {^12) and 
d^prodi^i'i) in Eqs. ( 11501 ) and (|151|) , so consequently for ^12 k 

'r{'di2) 1 (with leading particle). (153) 

i.e. the strong correlation has turned into a vanishing correlation. 

In case of running one finds the small angle behaviour again from the 
Born approximation. With p^p\'d) ~ {fi"^ /d) ln(ln(Pi?/A)/A) the direct terms 
without leading particle are given by 



<rL(^i2'^'0) ^ ^lnz(zln^-^ + l), ^ = lln^ (154) 

c?il(^i2,P,e) ^ ^rfJl(^i2,P,e), z' = lln^ (155) 

This yields for the normalized correlation 

, Inz' A/tlnln— , ^ 

r{^i2) ^ 1 ^ ^. 156 

m z V12 — K 
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One observes that the Bom term becomes less important at high energies 
in comparison to the fixed as case. After inclusion of the leading particle 
one finds again the dominance of the additional contribution with the result 

The correlation f{'di2) = P^^'*('^i2)/'"-^ vanishes for -012 — > k in any case 
either in first or second order of ln('(9i2//t). 



5.2.5 Comparison of various approximations for r(^9 
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In Fig. 7 we compare various analytic and numerical results on the (rescaled) 
normalized correlation r('(9i2). First there is the exact numerical result with 
and without leading particle at the finite energy for Y = InP/Qo = 5 ap- 
propriate for LEP (full curves). It shows a linear regime for small e < 0.3. 
The full result with leading particle vanishes for 'di2 ^ k = Qo/P, i.e. for 
e — > eo = ln(P6/(5o)/ ln(-PO/A) < 1, on the other hand without leading 
particle the correlation rises strongly in this kinematic limit. These results 
are obtained by solving the integral equation (|126|) numerically^ and with 



( |125| ). The inhomogeneous term is calculated using Eqs.( |12^ , p.29| ) with the 
exact p^^^'&u) from (^) where the fci-integral in (|129|) is derived also numer- 
ically. With increasing energy the cutoff eo — * 1 and the effect of the leading 
particle vanishes for F — *^ cx>. 



For these high energies the rescaled correlation \n r{i!} 12) / ^Jln{PQ / A) ap- 
proaches a limit for fixed e as given by Eq. ( |147| ). The linear approximation 
|/3e of this quantity (dashed line) represents the asymptotic result well for 
e < 0.6 but does not fully match the result for Y = 5. 



The nonleading terms in = l/{i3^Jln{PQ/A) given in ( |108|) provide a 
finite energy correction to r{i)i2). In leading order of e one obtains from 
143| , |108| ) the expansion 



lnr(^i2)/v/ln(Pe/A) = -/3e(l + -( + + ■ ■ •) (157) 

The first term corresponds to the leading power a/2 in ( p.37| ), the terms of 

^ Applying Simpson's rule to the integral in ( |126D on a grid of values for the logarithmic 
variables yi = In Ki, Xk = In^'fc in the range (Qo/'i?i2 < K < P;t9i2 < 5* < 0) we 
can express p^"^^ at each point (j/i,Xfc) by the respective values at the points (yi-iXfc), 
(yi^Xk-i) and (yi_i,Xfc_i). Starting with p*^^^ — S'^'^ at the lower limit of the integrals 
one can successively calculate p*^^^ for all points on the grid, finally for P, Q. 



38 



higher order in ( which amount to about 10% match better the exact resuhs 
at small e (upper straight line in figure). We conclude that the expansion 
of In p^'^^ in ( is rapidly converging for e < ^ and already the lowest order 
term which corresponds to a running with scale Q = PO/A gives a fair 
approximation for large opening angles. 



5.2.6 Correlation functions for quark jets 

In the study of e^e^ —>■ hadrons one is interested in comparison of data with 
predictions for quark jets. The correlation function Vai'&u) for arbitrary 
initial parton a can be derived from our results on gluons using relation 
Eq. ([l9|) . Including the leading particle according to Eq. (|125|) we obtain 



Pi'n^l2) = C„p(^H^12)-(Ca-a4L(^12)+2c,p«(^12) (158) 

= cld^^lM + 2cJ'\^^2) (159) 



and finally 



Ca<L(^?12)+2pW(^ 
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where the distributions on the r.h.s. refer to initial gluons. For sufficiently 
high energies p'^^^'&u) <^ d^Jodi^i^), then the leading particle effect becomes 
negligable and we obtain 

rai^l2) = l + -ifgi^l2)-l). (161) 

This equation means that the normalized cumulant correlation (r — 1) is 
rescaled by 1/ca (because F^^^ is multiplied by Ca and d^pJod by according to 
Eq. (p!9|) ) . Ultimately in exponential accuracy rq{'di2) Tgl-t^u) for infinite 
energies (a(P9) 0). 

In Fig. 8a we display as dashed curves the predictions on fa(i?i2) for quark 
and gluon jets from Eq.( |160|) where we insert the high energy approximations 



for gluon jets from Eqs. ( |142] , |140| , ^2D . The full curve represents the asymptotic 



limit for both types of jets. 

Next we consider the correlation function f with the alternative normal- 
ization fa = p'^^nl. This quantity is calculated from Eq.( |158|) using the 



39 



approximate formulae Eqs.( p.42]j4^j36| ) and Ua 



^a(e) = c„ ?/exp(-?/(2-^(e,2)) 



- (c:i-l)yexp(-2y(l- v^r^)) 

+ 2/5y%/(c„/(l-e)^)exp(-|/(2- v/r^)) (162) 



where y = 2/5y In (PO/A). In Fig. 8b we show the predictions for this quan- 
tity at finite energy for quark and gluon jets (ca = 4/9 or Ca = 1 resp.) as 
well as the asymptotic limit Eq. ( |145| ) . The finite energy effect at large e is 



smaller than for r('i9i2) where the leading term in the exponent is canceled. 
The shape of the asymptotic curve in the figure ~ 2/5(1 — ci;(e, 2)/2) is reached 
faster with increasing energy than its normalization. 



5.3 Fully differential connected angular correlations 
5.3.1 2-parton correlations for fixed 

This calculation shares common features of both cases considered earlier, 
namely the angular ordering, as in the single parton density, Eq. (|23|) , and 
integration of a product of the singularities as in Eq. (|26|) . We consider here 
the case of constant as ■ Results for running as will be discussed in the next 
section together with the ones for general order n. Integrating our resolvent 
representation, Eq. (|6lD over parton momenta gives 

{2)/n n ^ {2)/n n ^ i /" dK (KIk j, P QpK^ (2)/(-. x 
py{Qi,Q2)=9p'{^i,^2) + — / -^7^H-r^, )9k {^Kk„^Kk2), 

(163) 

with the resolvent I{x,y) defined in Eq.( [131| ) and the inhomogenous term 
given by the product of the single particle densities p^^^ (fi) = p^^^ (t?) / (27r'(?) 
in (|43D . Similarly to the previous cases, nested terms do not contribute to 
the leading behaviour in the high energy limit and will be neglected. If 
needed, they can be included without any difficulty. Saturating the angular 
integration by two poles 0A'fci(6xA:2) ~ we get for the connected correlation 
function, 

i^pfiau) f'" ^ {eK.AHe.,,)) + (i ^ 2)) . (i64) 
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As before tik was replaced by {fik^ in the 1^2 term) in all nonsingular 
expressions. In particular, QpK — ^ "^^i and QKk2 ~^ '^12- Note also that 
a{ki,k2,K) = ■di2 in this case since the minimal virtuality (emission angle 
Qpk) required for the parent K to emit ^2 is indeed controlled by if 
The upper limit of the Qrui integration is given by 'di2 and not by Opfci as 
one might have guessed from the angular ordering. This is a consequence 
of the arguments following Eq.(^) Q Similarly the lower bound of the K 
integration is controlled by di2 and not by Qpki- Integration similar to one 
performed explicitly in Appendix E gives 

vf{nr,n2) = 



- (2.)2.?.L W(2av//L0 + (1 ^ 2)(165) 

where I = In {^12/1^), Li = In {^1/^12) and ^1(^2) > ^12 respectively. Eq. (|165D 
is the new result for the full angular dependence of the two parton correlation 
function in DLA. Similarly to the density in the relative angle, Eq. (|136|) , we 
find a characteristic power behaviour at high energies (see Appendix G) 

which shows again the fractal (or intermittent) nature of the cascade. 

Again one can prove "backward compatibility" with earlier results. For 

example we show in Appendix H that integrating Eq. (|164|) over Qi,Q2 at 

fixed relative angle reproduces the second term of the resolvent representation 

for p^^^('i9i2), Eg . (|13CI|) as it should. Similarly one can integrate directly the 

asymptotic form, Eq.( [166| ), to obtain Eq. ( |136| ). Note that in both cases one 

(2) 

should include the disconnected term gp which is not negligible in the high 
energy limit. This has an interesting consequence: even though only the 
connected correlation F has a simple power behaviour, on the more inclusive 
level it is the density p'^'^\'di2) which is power behaved. 



^However the connected correlation F in Eq.(164) does not vanish only if Qpk^ (0pfc2) > 



'!?i2 hence the angular ordering along the cascade is preserved. 
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5.3.2 Running «s results for general order n 

Next we compute the cumulant correlation functions T^p\{^l}) fully differ- 
ential in the angles {Q} = {Qi, Qn) for running a^. In the DLA we require 
for the polar angle i)i of particle i and the relative angles "^im 

A,^im:^A/P. (167) 
We study the correlations again in the high energy limit with the quantities 

eL = H^i/^im)/ HP^,/A) (168) 

kept fixed within the allowed range (0 < < 1). We consider here the 
typical configurations where the various relative angles -^im between partons 
are of comparable order but with 

ICI « 1, C = ln(^.,M„)/ln(P^,^/A) (169) 

Note that this condition can always be satisfied for any e}^ for sufficiently 
high energies as 5ll^ = ln{'&ij/i^im)/ ((1 - ej^) \n{P'&i/A)). 

We start by deriving the integral equation for Fp'* from Eqs. 
First we carry out the momentum integral for Ap'' in (|2l|) , thereafter the 
angular integral is obtained applying the pole approximation: the integrand 
becomes singular whenever the direction of the intermediate jet K becomes 
parallel to any of the momenta ki because of the singularities l/i?i:fc. in 
^^g^K,prod- ^hc integral can then be written as a sum over n terms arising 
from the leading logarithmic sing ularities = ELi Apl with 

rP flK r^™' Hfpl^ 

(170) 

Each term has the singularity l/'d'i but depends otherwise only on the (^ — 1) 
relative angles VLk^k- between parton i and the others. This results from re- 
placing the slowly varying functions '&Kkj{QpK) by the value ^kikji^i) for 
The integral equation for Fp'' ( pOD is then solved -neglecting the con- 
tribution from the nonleading nested term- by the analogous decomposition 
Fp ^ = J2i=i Tp] with each term obeying the integral equation 

rK({.}) ^ Ag . 1 f £; f^a^m) [BL r£({^.})] • 

(171) 
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The bounds in the angular integrals of ( |1 70| ) , ( pTID correspond to the different 
regions T" and T' discussed in section 2: is taken as minimal angle between 
particle i and the others, "t^f^ is the maximal angle. For n=2, = -df^ = 
= di2. The lower limit in the K integration with angle df does not enter 



our high energy approximation. Furthermore in ( |1 71| ) the angular ordering 
condition requires for all relative angles dij < di. 

First we calculate aJ?^ for = 2 with d(2)(l,2) = p(i)(l)p(i)(2) and 
P^^\VL) from (|42|). In the integral over the first pole for ^Kki ~ we replace 
P^^\VL2) by ^^^-'(^12). The angular and momentum integrals can be integrated 



by part whereby terms of relative order 1/ \J\n{P^ 12/ K) are neglected (see 
( |167| )). Then we obtain 



where we have written a{P'di) = {P "d i) / a{P'd 12) ■ To derive Fp^ we note 
that the integral equation (|171|) is of the generic form (|28|) with h = "i^^F^^^ 
and pt = Kdi. Our inhomogeneous term (|172D can be written in the form 
Ap^ ~ exp(4/3-i/ln(P'(9i2/A)) just as in (|29D . So we can take the solution for 



h at high energies from Sect. 4. The prefactors a{Pdi) in ( p.72| ) as well as the 
"mixed scale" a{Kdi) in the integral equation (|171| ) influence the result only 
at the next to leading order of 1/ JlnPdi/ K (see Sect. 4). We obtain 



r?^(fii,fi2) ^ ^0|^e^P [2(3u{e\2,2)^HP^)) + (1 ^ 2). (173) 

where the prefactors ensure the correct limit F*^^-* — > A*^^^ for di — > = ■&12 
according to ( |171| ). With u{e, 2) ^ 2 — 3e/4 for small e (sufficiently large {^12) 
we obtain 

2(4^0^ fci \~^J +^ ^ ^ ^ ^ 

which has again an approximate power behaviour in the angles. In the a-limit 
( |37D of constant as one obtains back the previous result Eq. ( |166| ). 

These results can be generalized to an arbitrary number of particles. The 
inhomogenous term d^"'^ is now built up either from the product HiLi P^^^ (^) 
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or from combinations of correlation functions T^"^^ of lower order m < n. 
Although for rising n there is an increasing number of terms corresponding 
to these various different combinations the final result can be written in 
compact form as follows ( see Appendix I for more details) 



The intermediate scale Pi)i is defined by a{Pdi)2 = a^[Pdi)a'2 "^[Pd^) 
where is the minimal or because of ( |169| ),( p89D any other relative an- 
gle and is the maximal relative angle. The functions in ( |175D are 
homogeneous in the relative angles dki of degree p = —2{n — 1) and can be 
derived recursively. They are built from factors [kl] = l/i^li and read for 
n = 2,3 explicitly 

= F| = [l,2] (176) 
F,' = {3[12][13] + [23]([12] + [13])}/2; F| cyd. (177) 



In the linear approximation for uj{e, n) we obtain 




(178) 

with the by now well known power behaviour. Within our accuracy ( |167] ) 
this result approaches Ap^ for di — > df^ . 

The correlations implied by Eq. (|175|) are dominated by the coUinear sin- 
gularities from the gluon emissions. In each term i the particle i is connected 
with the initial parton direction (factor t?"^) and all other partons fc, / are 
connected among themselves (factor 'd'^^). This structure is illustrated for 
n = 3 in Fig. 9. The "star" connections in the left column originate from 
the first term in (^ or ( [^91| ) in App. I, the "snakes" in the other diagrams 
from the terms involving the connected parts. Our results for the QCD cas- 
cade correspond then to the phenomenological model by Van Hove ||2^ which 
builds up multiparticle cumulant correlations from products of two particle 
correlations as in Fig.9, whereas in the linked pair approximation ^9] the 
"stars" are absent (for a discussion, see pO|). 
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6 Local multiplicity moments 



6.1 Definition of observables 

In this chapter we discuss moments of multiphcity distributions in sidewise 
angular regions: (a) a cone of half opening angle 6 at angle i) with respect 
to the initial parton and (b) a ring of half opening 6, symmetric around the 
initial parton at angle i). The two cases correspond to dimensions D = 2 
and D = 1 of the phase space cells of the volume ~ 6^ (see Fig. 10). The 
unnormalized factorial and cumulant moments of the multiplicity distribution 
are defined as integrals of the respective correlation functions over the angular 
region 6) 

/(")(t9, 6)= f dh... dKp^^\h ...K) (179) 

C(") {d, 5)= I dh... dknT^""^ (fci . . . kn) (180) 

These quantities are conventionally normalized by powers of the average mul- 
tiplicity in the respective angular region 

F(^) = /('?)/n'?, C^'i) = 0^1)/^!^ (181) 

The factorial moments can then be calculated as the event average of multi- 
phcities in the region 7('(9, 5) by 

=< n(n - 1) . . . (n - g+ 1) > /n« (182) 

The derivation of analytic results for these moments from our correlation 
functions F*^") require some approximations. They can be avoided for mo- 
ments which are differential in one particle momentum. So we also consider 
the differential cumulant moments 

1 n „ 

h{k,S) = -y^ dk^...dknT^^\ki,...kn)5{k-ki) (183) 

where one particle is kept fixed at three-momentum k (or at angle d after 
integration over |A;|) and all others are counted in the region 7 around three- 
momentum k (or d) 0. We are particularly interested to investigate under 

^Such moments ("star integrals") have been also discussed recently by Eggers et al. 
13 
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which conditions the moments behave hke a power 



(184) 



Such behaviour is indicative of a fractal behaviour of the multiphcity fluctu- 
ations, also called intermittency. A jet with such a property can be assigned 
a fractal dimension 1331 , 103 



Dil- 



4>n 



n 



(185) 



For a random distribution of points in D dimensions one obtains a Poisson 
distribution in each subdivision, so for all 6 the factorial moments are 
1 and (pn = 0, so D„ = D. On the other hand, if in any subdivision with 
cell size 6 there is only one cell occupied and the others are empty, i.e. the 
distribution contracts to a point, then one obtains 0„ = n — 1 or D„ = 0. 



The general formula Eg. (185) covers also these two limiting cases. 



6.2 Integral equation for moments 
6.2.1 The equation 



One can derive QCD predictions for the moments (|179|) - (|183|) in two ways. 



In this Section we obtain and solve the corresponding integral equation. The 
asymptotic behaviour at larger angles can be also obtained from the direct 
integration of our results ( |16(j| , |178| ) for connected correlations. 

It follows from the previous Section that the connected correlation func- 
tions depend on a limited subset of all angles which characterize the con- 
figuration of the final partons. For example Tp^{fli,Q2) depends only on 
two variables -^i, ■^12. As a consequence the integral equation satisfied by 
moments themselves can be readily derived. 

We begin with n = 2. In this case integration (|183|) over the sideways 
cone 7('i9, 5) reads explicitly 

hP{n,6) = An f dd^2^i2T%{n,n^2). (186) 
Where both i = 1,2 give the same contribution. Integration of Eq. ([27|) over 
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the cone, Eq.( |186| ), gives for each i {i = 1,2),Q = {{}, 0), Qi = (pi), 

(187) 

Changing orders of the ■i?i2 and {K, -di) integrations produces two terms 

1 rP dK I ("^^ /n^ \ 

- / — / d,d{}ia\m 27r / ^9i2C?^?i2r ^^12) + 

i)^ JQo/6 K Jqo/K \ JQo/K J 

/ / ^4^^a\Kd)[2T, / ^i2rf^i2r?(a,r]i2) . (188) 

Inner integrals are nothing but the moments of the K jet Eq.( p,86|) . Adding 
the contributions from both i gives 

tf'^ jQo/5 K JQo/K 

IfP l\di,a\m)h^^\n^,5), (189) 

This result remains valid for arbitrary n as well. The essential property, 
namely emergence of the two terms in the integral Eg. ( |188| ), which can be 
reinterpreted according to the definition ( |186| ) is independent of the order of 
the moment. Equation ( |189| ) can be brought to the more symmetric form 
introducing the logarithmic density hp{Q,6) = 'd'^hp{Q,6). We get for arbi- 
trary n 

h^;\n,,6) = d^;\n,,5) + r ^ /' ^a^(ir^)4"\fi^,^) 

^I' dKr^d^ ,., 

jQo/S K JS %) 

— (n) 

Inhomogeneous parts dp are given by the corresponding integrals of, once 
iterated, disconnected contributions Ap ^ 

„q2 n, „ 

df{a,,5) = —Y. / dh...dK^f{ki,...K)5{k-k). (191) 
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In particular 

2 g 

dfin^,5) = 2nJ2 [ {&'Af\n,nu,p)) A2d^i2, (192) 
i=i •'''^ 

with Aj-^"* given by the i ~ th term of Eq.(^). Even though equation ( |190| ) 
has shghtly more comphcated structure due to the term h{^l^,ip), it can be 
solved in two steps - each step involving the type of equation encountered 
earlier. To see this note that at S = the last term vanishes and we get the 
simpler equation for the boundary value fi'p\ip) = h^p\VL^,'il)) , 

hP{5) = d^-\5) + r ^ /' ^a^K5)hPi^l.). (193) 

JQo/5 K Jkk W 

Equations of this type are satisfied by global quantities like the total average 
multiplicity n(PB), Eq.(^), or factorial moments in the whole cone (P, G), 
Eq. (|118|) . For constant as they can be easily solved by iteration 

d f ,1^,, . r . , . , .... .A , / P5\ 

X = m 



h^"\x) = — + a£ sinh [a{x - t)]S''\t)dt 



. 



(194) 

while for running as the asymptotic behaviour of the solution will be given 
below. Now equation (|190|) reads in two variables 

h'-;\n,6) = hP{6) + r ^ f'^a'{mh^;:\n^,6), (195) 

jQo/S K JS ip 

with the known boundary term. Hence the equation for moments, Eq.( |19CI| ), 
was reduced to, standard by now, integral equation with slightly more com- 
plicated inhomogenous term. 



6.2.2 The inhomogenous term 

The two step procedure leading to Eq.( |195D has a simple physical interpreta- 
tion. The inhomogenous term h^p\6) = h^p\Qs,S) is the cumulant moment 
in the maximal sideway cone which is allowed by the angular ordering, see 
Fig. 11. Therefore it should be related to the total cumulant moment over the 
full cone. In fact we shall prove that with x\ = ln{P6/A), Xx = ln(P'i9/A) 

c^^H^a) = dPiS) = ^^^^"H^a)|x.=.., (196) 
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where 

r(")(XA)= / A^^\n,...n^)dn,...n^, (197) 

Jr(p,i9) 

characterizes the full forward cone (P, Since the asymptotic behaviour of 
r is known, Eq. ( |196| ) gives immediately the corresponding expression for the 
inhomogenous term d. This determines h^^\5) via Eq.(193). Relation ( |196| ) 
explains also why h^'^\6) satisfies the integral equation which is characteristic 
of the global observables. To prove Eq.( p.96| ), rewrite the definition ( |191j ) 
using Eq.(pTD 

ZiT JQo/S li ^PK 

(198) 

Due to the angular ordering all internal angles in the cascade are negligible 
in comparison with Qpx, therefore Qp^ — Qpki = S = Qpcf\ , see Fig.ll, 
and integration over the initial parent dQ^ at fixed ki can be transformed 
into the integration over the final parton dQi^-^ at fixed K, see Appendix J 
for the detailed discussion of the n = 2 and n = 3 cases. We obtain 

= i In. ^^'mtPiS), (199) 
Zn jQo/s K 

where 

tPi6) = [ • • • , ^n)dn, ...n^, (200) 

is integrated over the complete forward cone of an intermediate parent parton 
{K, 5). On the other hand inserting Eq.(^) into Eq. ( |197| ) and interchanging 
orders of d^K and angular integrations one obtains 

.?)(tf)=r f r '-^aHmtfm (201) 

Comparing Eqs.( |199|) and ( |201| ) gives the relation (|196| ). 

Once the rule ( p.96| ) for the inhomogenous term d is established, we imme- 
diately see that the whole solution h^p\6) = h^p\\n{P6/Qo)) of the boundary 
equation ( |193| ) satisfies an analogous relation. To derive it, integrate Eq.( p.93| ) 
over S. One obtains 

^?(")(X,) = / Up\5)- = rW(X,) + / dy - dza\z)gn{z). (202) 



^The cone is defined around the direction of the first parton ki 
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Since r*^") is given by Eq.( |201| ) and t*-"-* is the source (the inhomogenous term) 



of the global factorial moment /„, cf. Eqs. (|117| . |T5D f\, we conclude that the 
solution of Eq. ( ^02| ) is simply 

g^''\Xx) = r^dy r dza\z)Uz), (203) 



A JX 



which can be also checked by the direct inspection. Therefore we have proved 
that 

h^''\xx) = ^ r dza\z)Uz). (204) 



27r Jx 

This relation determines immediately the high energy behaviour of the inho- 
mogenous term in the two variable master equation for moments, Eq.( |195| ). 
For, inserting Eqs.( |120| ) and (|36|) into ( p04| ) gives the following asymptotic 
result g(/i(")(XA) = hp (6)) 

(xa) = Vn exp (2n/? v^) , (205) 

and 

V„='^(^=] . (206) 

Inn ' - ' ' < ' 



2^^a{xx) 

with coefficients / and F„ given by Eqs. (|36|) and (|120|) respectively. 

6.3 The solution 

6.3.1 The leading contribution 



Integral equation ( |195|) is identical with the generic equation Eq. (p8D if we 
substitute pt K{}. Due to this universality, the asymptotic solution for 
running is readily obtained from the solution (|88D of Eq.(^). In particular 
the scaling variable is now 

e = ln(t9/5)/ln(P?9/A) (207) 

As was discussed in the previous Section, the asymptotic behaviour of the 
inhomogenous term is directly given by that of the global moments. With d 



Neglecting nested terms at high energy. 

^Equation (12C) is also valid for the running as with the replacement X X\. 



50 



of Eq.(^) given by Eqs. ( |2O50 and ( |206| ), and using the machinery of Sect. 4, 
we get the following remarkably simple result for the leading contribution to 
the unnormalized cumulant moments 

7^^"^ {6, ^, P) ~ exp (^2/3^1n {P^/A)uj{e, n)^ , (208) 

with the universal scaling function uj{e,n). This equation is one of the main 
results of this paper. It forms the basis of the following discussions and 
further approximations. The function uj{e) is determined algebraically by 
the saddle point condition, Eqs.( pO| - p^ ). A power expansion in e can be also 
easily obtained by the WKB technique proposed in Section 4.2, cf. Eq. (|107|) , 
as well as the 1/n expansion, Eq.( |ll(J| ), which approximates the full result at 
the percent level. 

The cumulant moments, Eq. (|180|) , are simply obtained by the integration 
of h = h/d"^ over the cone {D = 2) and over the ring {D = 1) for the two- 
and one-dimensional moments respectively. 

cW(^,5)~(^^j (^^j (209) 

After normalisation by the n-th power of the multiplicity in the respective 
angular region. 



one obtains 




n{^,6) = -^a{P^)f\-] ( — ) , (210) 



CW(^,5)~ ; ^ • (211) 



5 ) \A J 

where again D labels the dimensionality of the phase space cell. 



6.3.2 Nonleading terms 

The methods developed in Sect. 4 provide also the nonleading contribution 
to the asymptotic solution. In particular the normalisation of Eq.( ^08| ) is 
also available. Since however the double log equations by definition do not 
include all next-to-leading effects, the normalization of h constitutes only a 
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part of the nonleading terms. However we quote for the completness the full 
result given by the saddle point calculation 



h^'"\5,^,P) = CaVn-DnSxp (2/?^ln (P^9/A V(e, n)) , (212) 



with the prefactors Vn and Vn given by Eqs. (^) and (|206|) respectively. The 
factor Ca distinguishes between quark and gluon jets according to Eq.(|T9|). 
The factor T>n generalizes the n = 2 result, which can be also derived by 
the direct integration of the connected correlation function, see Appendix K. 
One obtains for the unnormalised cumulant moments 



cW(^,<5)=42-^vrc,P„P„(;) [— j , (213) 



and for moments normalized by the respective multiplicities 

^'"'(■M) ^ (i) (9 (^) . (214, 

We emphasize again that the double logarithmic approximation determines 
only the leading behaviour, Eq.( |208| ), and the normalization can be at best 
considered as the qualitative indication of the nonleading contributions. 

A second example of nonleading effects is the difference between the fac- 
torial and cumulant moments. Consider for simplicity the second moments. 
It follows easily from Eqs. (pi 3D and (|179|) that the difference /2 — C2 is non- 



leading 1^. Therefore, in our approximation, both factorial and cumulant 
moments have the same asymptotic behaviour and any distinction between 
the two can be discussed quantitatively only if the complete next-to-leading 
correction is known. 



6.4 Fractal behaviour of moments and Intermittency 

There are two angular regions where our solution reveals different behaviour. 

At small opening angle K/ P < 5 <^ d there are two large logarithms in 
the problem which can be chosen as \n.{d/5) and \n.{P{} / h). Both of them 

^°For large 5 {d/S fixed, P oo) it is of order a and for smaller 5 (e fixed, P — + oo), it 
is exponentially small relatively to the moments themselves. 
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are summed in Eq.( p08| ). At fixed e and P — > oo we predict the new type of 
scaling in tliis variable. However the effective exponent of the power law 




CW(79,5)~ - . (215) 



depends on 5 in a rather complicated way 

0„ = D{n - 1) - 2a{Pd){n - uj{e))/e, (216) 

since e depends on 6. The QCD cascade is not selfsimilar for small opening 
angles. 

On the other hand for larger opening angles with the ordering 1 <^ 
In {^/S) ^ In (P-d/A) we expand uj{e) and Vn{^) in powers of e (P-i?) In (■(?/ 6) 
and keeping the lowest nontrivial terms gives 



h^''\5,^,P)=V^r—\ (-) . (217) 



In the same approximation the cumulant moments read 

c(")(^,5) ^ 42"^7rPjif r (218) 



. A / 

where again D = 1, 2 corresponds to the one- and two-dimensional cells in 
the angular phase space respectively. We see that the normalized moments C 
have a power dependence on the opening angle 6. This result establishes the 
self-similarity of the QCD cascade when analysed in the angular bins of the 
size 5. Although the exponents in Eq.(|219|) depend on the polar angle i), they 



are independent of 6. Hence in this region normalized cumulant moments are 
simple powers of the opening angle. The well developed QCD cascade is 
strictly intermittent. 

Corresponding expressions for fixed are obtained by letting (P-i^/A)^'**-™-' 
— i> [Pd/Q^Y and a{P'd) — a as before. The result for the normalised mo- 
ment (|219|) translates automatically to the constant case as it does not 
depend on the scale Qq. it is infrared safe and free of narrow divergences. 
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We have derived here the power behaviour as the small e limit of the 
general result ( pi4| ). Alternatively these moments can be easily obtained by 
integration of the full correlation function as is demonstrated in Appendix K 
for n=2. We have also checked that the cumulant moment h^'^^ integrates to 
the global moment F^"^^ = 4/3 (see Appendix H). 



6.5 Discussion 

Only in the case of fixed the power law for moments is found in the full 
angular region at high energies (e fixed, P oo). The parton cascade is 
selfsimilar down to very small angles. 

Results for running ag are shown in Fig. 12 for the sidewise cone (D=2) 
and the ring (D=l) according to Eq.( pll|) . For sufficiently large 6 (small 
e) one observes the effective power behaviour in the opening angle, whereas 
at smaller angles (larger e) the moments reach a maximum value and bend 
downwards. This was explained in the previous Section. In Fig. 12c we show 
the energy dependence of the moments for D=2. At fixed e the dependence on 
energy P and angle enters through the anomalous dimension a{Pi)). Note 
that these calculations use the asymptotic form of uj{e) which vanishes for e — *■ 
1. A full inclusion of finite energy corrections would lead to an earlier drop 
and a vanishing of Inc*^"^ for 5 ~ Qo/-P or e ~ emax = ^^{Pd /Qq) / \ii{P§ / k). 
Numerical solution of the DLA equation for pi^^u) suggests that for P = 45 
GeV (Fig. 7) such corrections become important for e > 0.6. 

As was already mentioned, the exponent of the power dependence, at 
small e, does not depend on the opening angle 6. Even for running nor- 
malized moments show straightforward power behaviour in 5 albeit with the 
Pi) dependent anomalous dimensions. This result shows that the QCD cas- 
cade has the intermittency property considered by Bialas and Peschanski 
[0, although not directly in the limit 6^0 but for the opening angle 
much larger than the soft cut-off d{h./P'dy |^ that is for the well developed 
cascade. This cut-off tends to at high energy. Therefore only at infinite 
energies the cascade is intermittent down to to the very small opening an- 
gles. For larger e the formulae ( pil|) or ( pi3| ) should be used. The power 



behaviour at large opening angles 5 ~ ^9 (smalle) had been postulated before 
in a phenomenological model p| . 



.2, say) denotes the upper bound of the linear regime of oj{e). 
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The above power behaviour can be considered as resulting from the fractal 



structure of the cascade, i.e. the selfsimilarity of the branching process [17 



13] again only valid in the region where the running of can be neglected. 
One can associate a fractal dimension to the scaling property of the multi- 
plicity fluctuations and - using the above definition ( p.85| ) and neglecting the 
difference between factorial and cumulant moments - we find 

n 4- 1 

= -—^a{P^). (220) 
n 

This implies that the fractal dimension for all moments are determined by 
the anomalous dimension a{P'd) which governs the scale dependence of the 
average multiplicity, in particular D„ a{Pd) for large n The fractal 
dimension is independent of the geometrical dimension D of the problem, as 
a consequence the deviation from D is larger for D = 2 than for D = 1. 



7 Comparison with Monte Carlo Calculations 

In this section we compare our analytic results from the DLA with results 
derived from Monte Carlo methods which take into account fully the kine- 
matic constraints. We use here for illustration the HERWIG-MC program 



Iq) for the process e+e" — > uu with subsequent QCD evolution. This MC 
program is based on the coherent branching algorithm and takes into account 
the most important leading and subleading effects^. This allows for checks 
of the validity of our analytic results in DLA. Furthermore we can estimate 
the effect of cluster formation and hadronization on these results. 



7.1 Correlations in the relative angle 



12 



First we consider the correlation function r{e) normalized by the full mul- 
tiplicity in the forward cone. This quantity approaches at high energies a 

^^A similar result for the fractal dimension has been given before lia], but refers to the 



leading exponent of P-d/ A of the mmormalized moments in Eq.(21J), unlike our definition 



p5| ), see also |T|. 

^■^For a recent discussion, see also We take the parameters A = 0.15 GeV, mg = 
nig = 0.32 GeV. In our computations at the parton level we don't include the non- 
perturbative splitting of gluons into qq pairs at the end of the cascade. 
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limiting distribution in e, see Eg. ( P-45| ), 



^ ' 2/5(1 -^cu(£, 2)) (221) 



In^ 2 

where we dropped the contribution from the e-independent prefactor. As 
can be seen from Fig. 7.1 the results from the parton MC nicely confirm the 
prediction of scaling in e in the energy interval considered and in the range 
of e < 0.6. For larger e one observes a violation of scaling as expected from 
the influence of the cutoff Qq (see also the insert in the figure). The full 



curve represents the asymptotic limit of the correlation function ( p21| ). Here 
and in the following we adjust the absolute normalization of the curve to 
the data. Only the large leading terms in the exponent can be considered 
as reliably calculated in our present approach, the normalization is an effect 
of nonleading order and therefore may obtain large corrections. The dashed 
curve represents the prediction from the quark jet ( |162| ). If we used our 



calculated normalization the dashed curve would be shifted downwards by 
0.1 unit. 

Next we consider the same correlation in 79i2, but with the differential 
normalization of Eg. ( |123| )p^. We would again expect a scaling in e for this 
quantity, however, see Fig. 7.2, the data do not approach the scaling limit, 
Eq.(147). The finite energy corrections derived for quark jets from Eq. (|160|), 



see Fig. 5.3a, shift the peak of the distribution towards smaller e, as observed 
in the MC, but the energy dependence is predicted to be much weaker. The 
scaling violations in the MC are found most pronounced for large e, i.e. small 
relative angles 'di2- Apparently the disagreement comes from the normalizing 
function dpj^di'^i^) in Eq. (|123|) . This function is built from the 1-particle 
angular distribution p'^^^'d) which for small angles becomes very sensitive to 
the orientation of the jet axis. This will be discussed further below in the 
last subsection together with a possibly improved definition of this quantity. 



^"'in the MC the normahzing function d'^^J^^^e) (see after Eq.( |l22| )) has been calculated 
by "event mixing" , i.e. 'di2 is the angle between two particles in different events each one 
determined with respect to the sphericity axis. 
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7.2 Multiplicity moments 

In Fig. 7.3 we show the normahzed cumulant and factorial moments in 
dependence of e (by variation of 5 in (|89[)) for two primary energies and both 
dimensions D = 1,2[^. The theoretical curves for C^'^^ are computed from 
Eq. (|216|) . The description of the MC data improves with increasing energy. 



Again the normalization of the curves is adjusted. From our calculation in 
Sect. 6 we obtain for the D = 1 moment at e = 0, i.e. for the full forward cone 
of half opening 2^, the result C^^\0) = l/{3a{P^)ca) which yields 0.62 and 
0.78 at 45 and 900 GeV resp. for gluon jets (times 9/4 for quark jets). This 
is larger than the exact DLA result = ^(2) -1 = 1/3, see Eg. ([1191) . Our 
larger values at e = are due to the approximation 5 <^ -i? in the calculations 
for our typically small angular intervals. The strong decrease of data at small 
e is therefore not reproduced by our calculation at the lower energies (see also 
Fig. 6.2). Furthermore it is known that the DLA overestimates the value of 
the global F2 and that correct inclusion of energy-momentum conservation 



leads to a reduction of its value 35 



In the r.h.s. of Fig. 7.3 we show the factorial moments (F^^^ = C*^^^ + 1). 
At high energies C-^-* ~ F^^^ in the exponential accuracy of DLA. We have 
plotted our analytic results on C^^-* moments also with the data for the F^^^ 
moments. We observe that they fit actually better to the F-moments which 
indicates a peculiar cooperation of non- asymptotic effects. 

At large e the size of the D = 1 moments is underestimated. This may 
be related to our requirement ^^12 < S in our integration over the ring. 

The moments with D = 2 deviate strongly from the predicted shape at 
large e. The origin of this behaviour is presumably the influence of the Qq- 
cutoff at finite energies (see also the impact of an exact calculation for r{i)i2) 
at large e in Fig. 5.2), also a sensitivity to the choice of the jet axis (see 
below) could play a role - differently to the symmetric ID moments. The 
normalization calculated for e = 0, C(2)(o) = 4:/{3a{P^)ca) in Eq. ( pTl] ) 
would yield moments considerably larger than the MC-data. 

Next we test the universal behaviour of moments at high energies pre- 
dicted from Eq. (|216|) . In order to exhibit this universality we consider the 

In the MC we obtained the 2D data by selecting tracks from quadrangular areas with 
Ad = A(t> = 2(5 inside the ring. 
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quantity 



and in complete analogy the quantity F^^\ This quantity approaches for e 
fixed, P — > oo the scaling limit 

-C^") ^ 2/5(1 (223) 
^ 2/?(l-v/r^) (224) 

where the last expression is independent of n and follows in leading n ap- 
proximation for oj. In Fig. 7.4 we plot — C*(^) at different energies vs. e and 
the approach to the scaling limit Eq.( p^23| ) is nicely demonstrated. The 2D 
moments have about the same dependence on e at all energies although with 
an energy dependent normalization, but the dependence on e is steeper than 
expected from universality Eq.( |223| ). In Fig. 7.5 we show the MC results on 
F^'^\ The ID moments show a rather nice scaling behaviour in e with the 
expected e-dependence of Eq. ( |223| ) which sets in already at low energies. 

Finally we consider the n-dependence of moments. Results of analytic 
calculations were shown already in Fig. 6.3. In the quantity C*^"-* the lead- 
ing n dependence of uj{e,n) ~ n/(e) + 0(l/n) is removed. The residual 
n-dependence is shown by the curves in Fig. 7.6. They are reproduced 
reasonably well by the MC data. 

Remarkably, all angular correlations calculated here by the MC (except 
r{di2)) approach a scaling limit in e as predicted. Even very different quan- 
tities like f{e) and the D = 1 moments C*^^\ F*-^-* - they are calculated for 
particles in different regions of phase space, also f is a differential, C, F are 
integral quantities - approach the same asymptotic behaviour predicted by 
Eqs.( p21] , p23| ), disregarding the normalization. An experimental verification 



of this prediction would be an interesting confirmation of the universalities 
occuring in the QCD parton cascade. 



7.3 Hadronization 

So far we have compared our analytic calculations with the MC results at the 
parton level to check the quality of our approximations. Next we investigate 
how the process of hadron formation modifies the results obtained at the 
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parton level. For this purpose we apply the cluster model as incorporated in 



the HERWIG MC ||T6[. In fig. 7.7 we compare MC results at the parton and 
hadron level. There is good agreement in particular for small e (large angles) 
but some discrepancy towards the small angle cutoff. The agreement further 
improves with increasing energy. 

At sufficiently large angles the normalized quantities considered in this 
paper are infrared safe, i.e. become independent of the cutoff Qq. It appears 
that in this kinematic regime the hadronization effects are negligable also, 
i.e. the model realizes local parton hadron duality. This is an interesting 
new criterion for the applicability of this phenomenological rule. 

7.4 Limitations of the DLA 

One important problem is the choice of the jet axis which should correspond 
to the primary parton direction in the DLA. This we have chosen to be the 
sphericity axis (in some cases we also tried the thrust axis without appreciable 
difference). Some quantities, like the angular distribution at small angles 
with respect to the axis, are quite sensitive to the precise direction of the 
axis in an event, whereas other, integral quantities, like the multiplicity in the 
forward cone with half angle G around the jet axis, are less sensitive to this 
choice. This problem is illustrated in Fig. 7.8 which shows the distribution of 
partons p^^\t]) in pseudorapidity rj = — hi{'d/2) with respect to the sphericity 
axis. Also shown is the DLA result - using either the exact form ( |4T| ) or the 
high energy approximation (^2]) for p*^^)('i9) - which scales in x = Pd/A or 



\nx = ln(2P/A) —rj. We see that this scaling property is strongly violated by 
the MC data, in particular for large t] ^ rjmax or small angles ^9. Clearly, the 
jet axis defined in a multiparticle event does not coincide with the direction 
of the primary parton. Therefore the distribution p^^'^ ('d) is not well defined 
at small angles. 

This also explains the bad scaling properties of the correlation function 
r{'&i2), Eq.( p.23|) , observed in Fig. 7.2. Namely, the normalizing function 
is proportional to p''^''{'&i2)n{'&i2), see Eq. ( |128D , and therefore no 
good scaling properties can be expected for small {^u- On the other hand, the 
quantities f(e) and the D = 1 moments depend on the jet axis only through 
the large angles G or -d, furthermore there is a rotational symmetry, therefore 
these observables are stable against small fiuctuations of the jet axis. 



There is an interesting possibility to improve on this |21]. Instead of 
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defining an observable with respect to a global jet axis one calculates this 
observable with respect to the directions of all particles in an event in turn, 
each one weighted by the respective particle's energy. For the angular distri- 
bution, for example, one defines 

Pe\^) = I d%d'kj-^p^^\h, h)5{d - efc,fcj (225) 

This quantity is also shown in Fig. 7.8. As can be seen, the scale breaking 
is largely removed. In this paper we restrict ourselves to observables using 
the sphericity axis. It will be interesting to study the possible improvements 



using the modified definitions like (|225|) for our observables r('i?i2) and also 



the D = 2 moments to investigate the sensitivity to the jet axis. A further 
improvement would be the calculation of quantities like (|225| ) within DLA 
(see ref. for the azimuthal angle correlations). This would eliminate the 
global jet axis from both the theoretical and experimental analysis. 

Another serious effect of the DLA is the neglect of the recoil. The energy- 
recoil can be taken into account in the so-called MLLA [HI. Its effect for 



moments has been studied in ||2^ and is found to be typically of the order 
of 10%. There is not yet any scheme available to include the angular recoil. 
Its significance can be recognized from studying r('i9i2) in the full angular 
region < 79i2 < 2B. Note that in the DLA applied in this paper there is a 
peak near ^9i2 ~ but the correlation is minimal (r = 1) for < -^12 < 20 
because of angular ordering. On the other hand the MC calculation yields a 
second peak near — 20 (not shown) of the same strength which is due 
to the recoiling primary parton. Therefore the inclusion of the angular recoil 
in the analytic calculations seems to be the most important problem for the 
more quantitative analysis of angular correlations. 



8 Summary and concluding remarks 

Beginning from a single principle, i.e. the master equation for the generating 
functional of the QCD cascade, we have derived a variety of predictions 
for angular correlations between partons. In particular we have calculated 
the distribution of the relative angle p^^u), the fully differential two-body 
angular correlations r(f2i,f22) and subsequently also the general n— parton 
correlations. Multiplicity moments of general order in a restricted angular 
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phase space were also computed. All calculations were done for constant 
and running as , and the correspondence between these results was always 
demonstrated. The double logarithmic approximation applied here provides 
the correct high energy limit. Some nonleading effects, namely the ones 
which are included in the master equation were also investigated. 



This work summarizes and extends our earlier studies ||TT|-|jT^. Some 



of the results on multiplicity moments were obtained also by other groups 



21|, |2^, as mentioned in the introduction. 

An interesting aspect of the angular correlations considered here is their 
sensitivity to the soft gluon interference as the angular ordering prescription 
enters all calculations in an essential way. Indeed, in a recent study of polar 



angle correlations [^, using a quantity very similar to our p^'du), a large 
effect from the angular ordering has been noted - unlike in the case of the 
azimuthal angle or momentum correlations. 

Universal e - scaling. We have studied in detail the emergence of the 
power dependence on the scale characteristic in the given process. Whenever 
the observable in question depends on the single scale only, then, in the high 
energy limit, it is a simple power of that scale. This is the case for the total 
multiplicity and factorial moments in a forward cone, and for the angular 
distribution, e.g. n ~ Qi^'^^\ For fixed as (constant anomalous dimension 
a) one obtains a power law in Qi whereas in case of running the power 
law is violated by the scale dependence of a{Qi). 

For more complicated reactions, which in general depend on two or more 
independent scales, no straightforward power behaviour exists. We have 
found instead a scaling behaviour of a new type which applies universally to 
a large number of at first sight not related angular observables. In case of 
two scales Qi,Q2, the distribution of the relative angle and the second cu- 
mulant moment in the sideway cone, for example, show the universal scaling 

behaviour 

p~exp(^/5yh^/(e)), (226) 

where /(e) is a calculable universal function of the ratio of the logarithm of 
the two scales e = Up to some simple kinematic factors, the scaling 

function /(e) is the same for both quantities and analogous functions in 
e occur for the higher correlations. The general result, Eq. ( P26| ), implies 



that the simple power dependence on only one scale can be recovered in the 
"diagonal" scaling limit for fixed e, e.g. when Qi, Q2 oo,Qi = Qg- Indeed 
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then Eg. ( p26| ) reduces to 

p ~ Q^Wi)^(^) (227) 

with the anomalous dimensions depending on e. The same phenomenon was 
found for the single particle momentum distribution, Eq.(^), which also 
depends essentially on two scales. 

The "e-scaling", Eg. (|226|), has interesting experimental consequences. 



The observable quantity /(e) ~ In p/ a/Iu Qi should depend on the three 
relevant variables (for example 5, P) only through the single variable e, 



i.e. two variables are redundant. Preliminary results from LEP have 
verified recently this scaling prediction for the polar angle correlation [|1^ . 
The e-scaling holds approximately for the opening angles B in the range 
30° < e < 60°. 

Intermittency and fractal structure. Results of this work have a direct ap- 
plication to the intermittency studies which have recently attracted a lot of 
attention. This property of the parton cascade is obtained from the general 
formula ( p26| ) in the limit of small e, i.e. sufficiently large opening angles, in 
the linear approximation of /(e). The suitably normalized observables then 
behave like r ~ Q2"'^^^\ i.e. there is a power law in one scale Q2 (say i!^/6) if 
the other scale Qi (say P'd) is kept fixed. So the perturbative cascade is in- 
termittent in the appropriate angular variables and intermittency exponents 
are available analytically, c.f. Eq. (|216D . For fixed as the power behaviour for 
the angular correlations is actually approached in the full angular region for 
high energies. 

This power dependence can be related to the selfsimilar structure of the 
partonic branching process. One can assign a fractal dimension Dn to the 
normalized multiplicity moments and this is found to be given in terms of the 
anomalous dimension of total multiplicity a = 70 = ^JqoJtx, see Eq. ( p20| ). 

Monte Carlo tests. We have also confronted the analytical predictions 
with Monte Carlo simulations which don't have many of the limitations of 
the double logarithmic approximation. The results are rather satisfactory. 



One indeed sees the approximate e-scaling (p26D and the shape of the scaling 
function is approached with increasing energy. The factorial moments are 
closer to the asymptotic limit than the cumulant moments. In the small 
e region the correlation functions and moments are independent of the pt 
cutoff Qq. In this region we find the hadronization effects negligable in the 
Monte Carlo (HERWIG) calculations. On the other hand deviations occur 
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for particles with small angular separation. 

To conclude, there exists a rich spectrum of analytic predictions from the 
perturbative QCD for multiparton correlations. The double logarithmic ap- 
proximation provides the leading order predictions applicable quantitatively 
at very high energies (~ 1 TeV) but some interesting experimental tests of 
the theory in its relatively unexplored regime are suggested already at present 
energies. 
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Appendix A 

In this Appendix we shall derive asymptotic form of the solution Eq. (|86|) . 
Consider first the case of constant as ■ This will allow to rederive the results 
obtained also directly from the correlation function, cf. Eq.( p04| ), and will 
provide some orientation about the relevant range of variables contributing 
to the integrals in question. Inserting moments, Eq.(0), into Eq.(^) gives 
for the leading term 



with the source term dn{K,6) = VniKS/Qo)^"". After changing the variable 
in the inner integration to u = In we can extend the upper limit to the 

infinity. If the contour 7 is properly chosen this does not affect the leading 
behaviour of the result. One obtains 




(228) 




(229) 
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P dK 



exp 



1 p 

- In — (Vs2 + 4a2 - s) 

2 K 



K5' 



IQo/5 K 

Similar procedure for the K integration gives the pole in s 



ds 



iTii + 4a2 + 4s2 - 2s) - na 



-(230) 



exp 



^ln|(Vs2 + 4a2 + s) 
2 



This pole saturates the contour integral giving finally 



(231) 



where n = Qq/ P . Normalization T>n follows from Eq. (|204| ). We obtain 

V 



7r2"+in 



(232) 



with F„ = lim,^oo/n(a:)/ntoi(x)". With = 4/3, cf. (|TT9|), Eq.(^) 
reproduces Eq.(|304D for n = 2. 



Runnung as ■ Similarly to the previous case the saddle point s* in Eq.(|8^) 
scales as the coupling constant hence s* ~ l/vmP. Therefore it is sufficient 
to use the appropriate asymptotic forms of the confiuent hypergeometric 
functions. They can be derived from the well known integral representations 
of the solutions of the confiunt hypergeometric equation [Q. We obtain 



lim -^(o-, b, z) 

a— >+oo,z^+oo,6=co?ist 



1 



exp 



^ 4 
■A 



exp 



Y)^(^' + ^A)-^ (233) 
In 



+ 1) 



A VA 



lim $(a, 6, z) 

a—*oo,z^oo,b=const 



exp 



z{z + A)-z) 



exp 



(234) 



A, + + A 

— m 

2 2 
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where A = 4a — 2b. This imphes the following asymptotics for the A and 
functions required in Eq.( 



jF(s, x) ^ 

Ais,X) ^ 
9±{s,x) 



2(3 



X 



27r/32 l+px g--|3_(s,A) 



s2 + 4/5Vx ± s, 



(235) 

(236) 
(237) 



where = y 1 + 4^. Now Eq.(p7D can be rewritten as 

K[d,i),P) = j exp(sln-)J^(s,ln— )M„(s)expw„(s,ln— ) — , 



where 



PS. 



Un{s) exp Wn{s, lu ^) = l!nQo/A '^^ ^) ^^P ^"(■5, CT 

2/32 



A 

vJs, a) 



(238) 



(239) 



a 



In (— ^+(s, a)) + 2nPy/a, 
s 2p 



l+Po 



and the prefactor resulting from the inhomogenous term (i„(cr) 
follows from Eqs.(Mll20|) and (|3|). 



(240) 



Vn exp {2n(3Ja) 



a{a)Fn ( f 



2™ \2ja((r 



Saddle point condition for the a integration daVn{cT, s) 

2np 



(241) 



gives 



which has the solution 



2\2 



(242) 



(243) 
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At the saddle 



Vn(s, a*) = —(1 — n^) H — — Inn = Wn(s) 
s s 



(244) 



which is independent of the P6 and corresponds to the insensitivity of the 
leading contribution to the integral ( P39| ) to the upper limit of integration. 
The prefactor reads 
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Ur,\S 



-u„\a 



, ^ 1 — - 



(245) 



We now follow similar procedure for the s integral. Writing 

/i„(5,^9,P)~ / f/„(s)expW^„(s,L,/)^, 
we have from (|3|,|3|,|3|j24|) L = In (P^/A), I = In {PS /A) 



Wnis,LJ) = siL-l) + 

ho (s.L^ ^ 



lLg_{s,L) + ^ln[^g_,{s,L)]-f^{n'-l), 



2(3 



and 



Un{s,L) 



u„\s . 



(246) 

(247) 
(248) 

(249) 



Saddle point condition, dsWn{s, LJ) = 0, gives after some algebra Eq. (|9^) 
with z = s\/L/I3 . In addition Eq.( p48| ) can be simplified at the saddle point 
reducing finally to 

Wn{s\ L, I) = 2pLu{e, n), (250) 



and the overall prefactor 



1 



27iw:(s''' 



(251) 



which after some algebra gives Eq.(^) for Vn ■ It is easy to check that in 
the constant limit our Eq.(p8D reduces to the result ( p31| ) derived by the 
exact integration of the resolvent representation (|229|) . 
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Appendix B 



In this Appendix we shall derive the constant as limit of the energy moments, 
Eq. (|49|) . This is done by leting /?, A — oo while keeping the ratio P^/X = a 
constant. In this process the first parameter and the independent variable of 
the hypergeometric confluent functions become large. Appropriate asymp- 
totic forms of these functions are derived in Appendix A. Using Eqs. ( P34] ) 



and ( P35| ) one obtains for the ingredients of the first part of Eq.(49) 

$(a - 1, 1, sX) ~ {sX)-^/\sX + 4(3^3 - 2)-i/^exp {(3^3 - 1/2) 

)-3X^ 



exp 



exp 



^3X{3X + 4(3^3 



- 1 In 



sA + JsX + 4/3^3 - 2 



F{a,2,3{X + X)) 



4/32 



-1/4 



p{X + X)V2 



exp 



3{X + X) 



4/32 



\ 3^{X + X) 



1 



3^{X + A) 

2/32 

exp 



((3\ 3\x + xy 



In 1 



4/32 



\ s2(X + A) 



All symbols are as in Eq.p9|) except that n was replaced by the continuous 
variable 3. Constant Og limit is performed at fixed values of the kinematical 
variables. In particular, X/X << 1 and sA >> 1 in this limit. Expanding 
above expressions and keeping track of the 0(A) and 0(1) terms in the 
exponents we obtain 



$(a - 1, 1, sA) 



A 



3 + Vs2 + 4a2 ^ 2 



2V sVs2 + 4a2 



exp {f3'/3 - 1/2) 



(252) 



exp — — (Vs2 + 4a2 — 3) exp 



2Vs2 + 4a2 



exp 



F{a,2,3{X + X)) 



2p\ s + Vs2 + 4a2\ 

In z exp 



exp [ — — In (sA) 



2a2 



2s 



/3(X + A)V27rVs2 + 4a2 



Vs2 + 4a2(s + Vs2 + 4a2), 
/ a^X \ 

: exp 



(253) 
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exp I — - — (Vs^ + 4a"^ + s) 1 exp I — ^ I exp I 
exp In exp 



2s / V sVs2 + 4a2(s + + 4a2) / 

Where we have aheady replaced the 0(1) ratio by a. Finally, insert- 
ing these expressions into Eq. (^91) and using Stirling formula for the gamma 
function, we see that the divergent, leading terms canel and the first non- 
leading corrections give the finite result which reproduces the first part of 
Eg. (|5^) ■ Second part is obtained in the analogous way from the second term 
of Eq.(ll). 



Appendix C 

We shall discuss here some details of the derivation of the simplified resolvent 
representation, Eq.(|86D. First we will derive Eg . (^5]) . To this end recall the 
resolvent solution for the two dimensional inclusive density, Eq. ([77|) . Per- 
forming the steps outlined there one obtains with the running ctg 

Formula for the momentum distribution follows from Eq. ( |254D upon the an- 
gular integration. Changing orders of the and \1/ integrations, recalling the 
definition (following Eq.([79D) of R and using the substitution rules (|^) gives 
Eq.(p5|). In the reduced variables x = In |, t = In 2; = In |-, r = In A = 

In ^ it reads 

hn{x,t,\) = dn{t + \)+ f dz f dTR{x-z,t-T,z + T + \)drdn{T + \). (255) 

Jo Jo 

Integrating by parts and noting that d is independent of z we obtain 

hn{x,t,X) = dn{t + X) + f dTR{x,t - T,T + X)dn{r + X). (256) 

Jo 

where we have introduced 

R{x, t-T,T + X) = r dtR{x - z,t -T,z + T + X). (257) 
Jo 
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Because R has the Laplace representation, Eg. (|87|) , which is dominated by 
s ~ 0(1/ In one can set to the leading order diR{xi,X2, A) = d2R{xi, X2, A). 
Therefore finally 

R{x,t-T,T + X) = [ dzd^R{x-z,t-T,z + T + X) = R{x,t-T,T + X), (258) 
Jo 

since the contribution from the upper limit z = x is negligible This 
concludes the proof of Eq. (^) which can be also simply summarized as the 
dominance of the lower limit in the integration in Eq. (^) . 

Appendix D 



In this Appendix we shall derive the asymptotic behaviour, Eq.( |120|) , of the 



total factorial moments, and the recursive relation ( |121| ) for constant . 
First, we simplify the resolvent solution for /„(X). Changing the integration 
variables Z and ( in Eq. (|70|) to t = Z + 3 and v = Z + ( and interchanging 
the order of t and v integrals one can perform the t integral to obtain 

/„(X) = bniX) + a j"" sinh (a(X - v))hn{v)dv, (259) 

with X = y + S = In (PQ/Qo). We now prove Eq. (|120|) by induction. For 

n = 2, 

b2{X) =< n B2 exp (2aX)/2l (260) 

Assume now that 

/I \ 

6„(X) = 5„ -e'^^ . (261) 



.2 

Inserting ( |261[ ) into ( |259D one easily verifies ( |12(j| ) with 

n 



ft 

Bn. (262) 



— 1 

Since &„(X) is built recursively from the products of moments of /ower order 
(for example &3(X) =< n >^ +3 < n > (/2(X) — 62(-^))) with the same 

^^We were also consistently neglecting terms proportional to d\R{xiTX2,X) since they 
are nonleading comparing to the derivatives with respect to the first two arguments. 
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combined anomalous dimension, the form ( p61| ) follows. This completes the 
inductive proof of Eq.( |12(]| ). 

In order to derive now the recurrence relation ( |121| ) we rewrite the master 
equation, Eq.(|^), for the generating function of the global moments, i.e. we 
set u{K) = u. In the high energy limit the nested terms are negligible 
which amounts to puting u = 1 in the exponent of Eq.(|^). With these 
simplifications Eq.(|^) reads 

Zp{u) = exp{J d^KMp{K){ZK{u) - 1)). (263) 

Differentiating n times one obtains 

""1 /n - l\ 

f4X) = J2\^ ^ yk{X)g^n-k{X)+g^{X), (264) 

where we have introduced the cumulant moments 

g„.{X) = J (fKMp{K)Zp\u)\u=i = f^{X) - h^{X). (265) 

Using the last equality and Eqs (|120| , |261|) together with Eq.( p62|) to eliminate 
Qm in favour of gives the recursion Eq. (|121|) . 



Appendix E 

In this Appendix we shall derive the series representation Eq.(|135D for the 



distribution of the relative angle §12 for constant . The contributions from 

the product and nested terms will be calculated separately. 

A. Product term. Inserting Eqs. (|131|) and ( |132|) into Eq.( |13(]| ) and expanding 



in powers of a one obtains upon integration over \E' {k = Qq/P) 

Iprod = Pprldi^u) - 4?od(^12) (266) 

2a2 ^ a'^+' 




with 



^ _^ ^^^^^ 



°^ 4"^ _ 1 ry 
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Integration over z is proportional to the Euler beta function 

z/) = /q^ t^~^{\ — tY~^. Changing the orders of the r and m summations 
we obtain 

Where y = y'^/4: and y and z are defined below Eg. ( P-35| ). The second term 
in the curly bracket gives after summation the inhomogenous term dprod, 
therefore the inclusive density itself is given by the first term only. 
B. Nested term. Calculation of this contribution is entirely analogous to the 
previous one. After expanding dnest, Eq.( |133|) and the resolvent ( |131| ) one is 



left with elementary integrals. Final result for the nested part reads 



Ai^u) = ^(f: T^'''h^+i{z) - sinh {aHA2/K))] . (269) 



Adding both contributions gives the result ( p.35| ). 

Appendix F 

In this Appendix we integrate the formula, Eq.( |135| ) to obtain the total 
factorial moment /2 as given by Eg. ( |119| ). Integration of the second term in 
is trivial hence we write 



^ d^i2P^^^ {^12, P,Q) = J (a In (^-^^ ^ - 2 (^cosh (^In -^^ - 1^ . (270) 

with J{aX) ,X = In (G/k) defined as the integral of the first part of Eq. ( |135| ) . 
Expanding modified Bessel functions one obtains 

J{aX) = VV-l_-4 _/ (271) 

°° ( n \^'\2m+2fc+2 

i'o to(2^ + 2A: + 2)!' 
which can be resummed 

oo y2Z+2 I 

AY) = E 7^7-^ E S^'"^' (272) 

1=0 l^*- '^J- m=0 

= ^(cosh(2F) -cosh(F)). 
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Inserting ( p73| ) into ( |27CI| ) gives 

r P^^\^i2)d^i2 = \ (cosh(aX) - 1)' , (273) 
in agreement with Eq.( |119| ). 

Appendix G 

We derive here the asymptotic form, Eg. ( |136| ), of the correlation function 
at high energies k = Qq/ P ^ 0. In this limit 



y = 2 
with 



ln(^i2/«:) 



'12y 



oo, z = 2a Jin {^u/n) In (Q/^u) oo, (274) 



\\ in(eM 

z/y = a In (G/i^is) = = const. (275) 

Using the sum rule, 

oo .2 

J2y'h+,{z) = z-''ey^/' e-^^'/(2^)4_i(r)r-'^cir, > 0, (276) 
one can rewrite the sum in Eg. ( |135| ) as 

p(2)(^i2,P,e) = re--'/^'lo{r)rdT + e-^y/' f e^' '^^ h{r)rdr) 

O^V\2 \ Jo Jo / 

(277) 

with the first term 0{{'di2/ n)'^"') dominating the high energy behaviour. We 
have neglected the terms 0(('i?i2 //«)")• Since the integrand in the leading 
term is well behaved at large r, we can extract the large P limit by setting 
the upper bound to infinity. This gives 

p(2)(^9i2, P,Q) = ^ (—] r e-^' Ioiau)udu, (278) 

\ K J Jo 

which can be integrated exactly to give finally 

p(^)(^,2,P,e) = -^f^)"e-V4. (279) 
2'i9i2 V ^ / 

This is identical with Eg. ( |136| ) after taking into account Eg.( |275| ). 
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Appendix H 



In this Appendix we shall perform several consistency tests of our results for 
the two-parton correlations. 

First, we show that the form, Eq. (|164|) , of the fully differential connected 
angular correlation is consistent with the Eg . (11301) for the correlation in the 
relative angle. It suffices to show this compatibility for the connected parts 
only, since the relation between the direct terms was already established (cf. 
Eqs (^) and ([132]) ). To this end we define 

r^'H^u) = I dn,dn^5{d,2 - Qk,k,)T^^\^u^2). (280) 



with V'^^\ni,Vl2) given by Eq.( |164| ). Due to the symmetry between partons 



if is sufficient to analyse only the first term (Fi, say) in Eq.( |164| ). Choosing 
the polar coordinates of the dVL2 integration relative to the direction of the 
momentum of the first parton one obtains 

F(2)(^i2) = ^Jd^iJ &k,k,dQk,k,d(j)ur?\ni,ni2)S{^i2 - 0^,,,), (281) 



which gives after performing integrations over &Kki (in ( |164| )), 6„j„2 0i2 



2a^ /■© dK 

sinh(aln(irt9i2/Qo)) T /(^, 



2n^ /■© dK 

T^'\^i2) = ^ — (cosh(Gln(i^^i2/go))-l) (282) 



i}i2 K 6 



pfci 



This is identical with the resolvent solution ( 130 ) after changing the integra- 
tion variable Qpk^ to \E' = 0'i9i2/0pfci- 

As a second consistency check we integrate the cumulant correlation 
F*^^^(fii, VL2) up to the multiplicity moment /2 for the full forward cone which 
is known [| to be /z/n^ = 4/3, cf. Eq.(lT21|), with either p^'^\di2) or 
h^'^\VL,5) and S-'^\d,5) as intermediate results. 

In the first case we derive the cumulant correlation in the relative angle 
by integrating Eq.( |166|) . One obtains 



F(2)(^12) = j dn,dn25{^,2-Qk,M)^^^\^l^^2 
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a re 2nd^^^_^ ^^ ^ . + (1^2) 



Here we used again the pole approximation with the same results for both 
terms. With the results Eqs.( P83| , |136| ) for p^'^\di2) and d!p}od{'^i'i) we recover 



the relation ® V^^X^^^) = P^^\^i2) - <L(^i2). 

We can now calculate the multiplicity moment /2(5, 0) in a cone of half 
opening 5 averaged over the forward cone B and the corresponding normali- 
sation n^((5, 9) 

= TrfWH^is) ^ )^ (284) 

Jk 6 k 



n\6, 0) = / rf^i2CL(^i2) ^ 7(-)'^ (285) 



' ■^2<K^l2)^^(^)^^ 

where the approximation is for 6 ^ k. For the normalised moment we then 
obtain F2{5,Q) ~ |(f )^ • This reproduces F2 = | in the full forward cone 
(5 — s> 0). So we find here that the density p^'^\'&i2) and not the cumulant 
correlation T^'^\'di2) follows a power law, although the difference between 
both quantities disappears for 'di2 -C 0. 

In the second case we note that the cumulant moments h^'^\il,5) and 
c(2)(fi,5) derived from T^'^\ni,n2) are power behaved, cf. Eqs.( p0^j306| ), 
and not the factorial moments. To check the consistency with the previous 
results we calculate in analogy to Eqs.( p84] , ^85| ) the cumulant moment in 
cone 6 averaged over the full cone 0. 

c(2) (5, 0) = r r dnh^^^ (fi, (5) , (286) 

The first term appears because of the requirement 6 < i). Using our result 
(|304D for h^^\n, 5) we obtain for 5 > k 



c(^)(5,0)^i(V[|(|)t-l]. (287) 

With the normalisation ( p85| ) and the relations (^) we find again the consis- 
tency with the previous result on f^'^\S, 0) in ( p84|) , and therefore the correct 
integration to the full cone with F2 = |. 
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Appendix I 

Here we discuss the derivation of Eq.( |175| ) for Fp^ for general order n. First 
we have to calculate A^") from ( |170|) . In case the inhomogeneous term S"^"* is 
built up from the products of p*^^^ (i) we apply again the pole approximation to 
the angular integral which yields n terms, after partial integration in leading 



order of l/-^ln(i\:^9f /A) for z = 1 

. (n) _ /^/" [ 



(288) 



In the exponent and in the denominator we expand 

1 



In [Kdu/A) ^ {mf / + ^-5^^ + . . .) (289) 
and neglect the 5-terms because of (|169|) . Then, after partial integration, 



we obtain in leading order of l/\JP'd'i-/K and with the corresponding terms 
i > 1 added 

where i^" = l/(^fi ... ei/?+M ••• O- 

We move on to the case with S-"^"* as products of F^"^). Consider first S^^ = 
p^^^(l)F'^^)(2, 3). In the derivation of Ap\ in pole approximation we replace 
u}{e\2,2) ^ u}{-6li,2) ^ 2 in F(2)(2,3) of Eq. ([T73D , the same replacement 

(3) (3) 

applies in leading approximation of the angular integral for Ap2 and Apg 
so that again the form in Eq. (|290|) is obtained, but with different F". The 
complete inhomogeneous term from (^ has the structure (without couplings 
and exponential factors) 



d^""^ '-:^2:^ + :^n\:^ + iK^]+''y''^- (29i) 




Again, the angular and momentum integrals yield a result as in (|290|) , except 
for different functions F^, see Eq.( [177| ). In the same way the result ( |29(]| ) 
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follows for an arbitrary number n of particles and F" is again homogeneous 
of degree p = —2{n — 1). With the inhomogeneous term Ap] from Eq. (|290|) 



which is of the form (p9D we can then solve the integral equation (|171|) for 
Fp] using the results of section 4 which yields finally Eq.( p.75| ). 



Appendix J 

In this Appendix we shall discuss assumptions leading to Eq. ( |199|) . For n = 2 
Eq.([l9TD reads 

S^\6) = J d^Kdn2Mp{K)pi\nKc)pi\^2). (292) 



According to the definition, Eq . (11831) , we have chosen the cone axis C to 



coincide with the direction of the first parton ki. Integration over Q2 gives 
dP{6) = J d'KMp{K)p^^\QKc)n{KQKp). (293) 

Remaining integral is dominated by the singularities in Qkc due to the an- 
gular ordering {Qkp > ©Kfci = Qkc)- Using the pole dominance {Qrc ~ 
0, Oi^P ~ Bpc = f^) allows to integrate over dQx = dQxc and we obtain 
Eq.(|I^ since t'-^\K5) = n'^{K5). 



The case of n = 3 is more instructive since it illustrates other complica- 
tions due to the multiparton kinematics. We shall discuss it in more detail. 
The subsequent generalization to arbitrary n will be straightforward. First, 
the product term d'^^J^^ is given by Eq.(|D, and consequently one obtains 
Eq.(pUTD with 

tg)(^) = n^K^) + 3n{K^)g2{K^), (294) 



where 



g2{K^)= f Fg^(fii,fi2)rfl^irf^^2. (295) 

JiK,<f) 

We move now to Eq. ( P-99| ) for d^^\ Contribution from the product of the 
single densities to rfp^id can be treated in the full analogy to the n = 2 case. 
One obtains Eq.( |199D with the piece of t^^\ Similarly calculation of the 
contribution from the p^^-'(l)F*^^)(2, 3) term is simple. One obtains 



d'KMp{K)p^^\QKc)g2{KQKp), (296) 
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which due to the dominance of the Qrc ~ region gives Eg. ( |199| ) with a 
third of the whole mixed contribution 3ng2 to t^^\ 

The last two contributions are equal up to the relabeling 2^3. Consider 
the p*^"^''(2)r'^^)(l, 3) term in dprod- Its contribution to is 



d^KMpiK) ( c/fi2c/^^3pS^^(2)rg^(l,3). 

7-7 



Integration over ^2 gives 



d^KMp{K)n{KQKp) 



(297) 



(298) 



There are two regions which dominate the remainimg integrals: A) K\ \ki{C), 
this corresponds to the sequence P K ^ ^{C) 3(0p/^ > Qrc > 
Qcks) and B) KWk^ corresponding to the process P K ^ 3 ^ ^{C) 
with the ordering Qpx > ^Kk^ > ©cfcs- In each of them the connected 
correlation function depends on different set of variables. However in both 
cases dQxd^s integration is equivalent to the complete integration over the 
phase space of the final partons, which reproduces the total cumulant g2. To 
show this let us recall the generic expression for g2, Eq .( p95|) , if the underlying 
process (S) has a two-step structure K ^ ki ^ k2. Then r^''(f2i, ^^2) = 



r^Hfi/rfei,fii2),0itfci > ©12) and Eg . (12951) becomes 



92 



dVL 



Kk, 



dn,2r^> 



12 



(299) 



where k = Qo/K and limits of the integrations apply to the corresponding 
polar angles. 

Now, in the case A) we obtain the following contribution to Eg.( |199| ) 



(A) 



P dK a" {KQkp) dQK dQ^ , , 



Qo/S K 



'^^ic'^'ckS^i'^KC-, Qcka] 



fi2 f)2 



(300) 



dK 
271 Jqo/s K 



a\K6)n{K6) 



^dQK 



c 



e 



277 



KC 







e 



277 



C3 



^KC^ 



ck'i'- {Qkc, 
P dK 



e 
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which gives the half of the ng2 contribution to t^^'^ in the virtue of Eg. ( |299| ) . 
In the case B) corresponding contribution to reads 

1 rP dK r 
i^) = ^ -rra\K6)n{KS) / driKdnsV'' {QKk,, Qck,) (301) 
Ztt JQo/s K J 

which can be rewritten using the pole approximation for the inner vertex 
1{C) -^3:eck,--0,QKks-QKC 

1 dK r^Kc 

(^) = TTT^ / -rra\K6)n{K6) / deKc27T / dec327cr'' [QkcQc 

ZTTO^ jQn/S K Jk Jk 



(302) 

Comparison with Eg. ( p95| ) shows that this provides the second half of the 
ng2 term to t*-^^ Adding the contribution from the p'-^''(3)r'-^''(l, 2) gives 
finally complete mixed term 'ing2- This establishes Eg.( |199| ) and hence the 
derivative rule (|196D . 



Appendix K 

Here we derive the second multiplicity moment by the direct integration of the 
cumulant correlation function. We first calculate the differential cumulant 
moment h{VL, 5) with one particle kept at angle Vl and the others within a 
cone of half opening 5. We start again with = 2, first for fixed a^, using 
Eg.([l66D 



with fii2 = (5i2,0i2)- In the second term we use -(92 = "&\ + '(}\2 + '^i'&i2Cos(t)i2 
and apply the pole dominance approximation for the i9i2 integral, i.e. we 
set di2 = in the nonsingular factors. This returns ^2 = and makes the 
second integral egual to the first one, so we obtain 



/ dQ2 r(2)(fii,fi2 



. a , 

2^' 



^ §i2d^i2d<Pi2 d 



,92 

'^12 



2a 



nn2 



(303) 
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The cumulant moment, calculated from the definition ( |18CI| ) has a compli- 
cated boundary for the fli,Q2 integration. With the approximation 



c(2)~ J ^^^h^^\n,6)dnc^7c6^h^^\n,s), (305) 

we obtain 

The average multiplicity in the cone 7(19, 6) can be derived within the same 
approximation from Eg . (^31) as 

^(^,^) = ^(^)^(^r, (307) 



and we obtain for the normalised moment from Eg. ( |181| ) 

C(2)(^,5) = f (308) 
3a 0^ 

In the case of running ag one obtains in complete analogy by integrating Egs. 

^^\^,6) ^TT6^h^^\n,6), (310) 



1741 H) 



c 



4 ,^92 



The consistency of these results with the case of fixed as is easily estab- 
lished by letting (3,X ^ 00, = a, (see Eg. (|37|) ). 

We have checked the guality of our approximations in case of fixed a^. 
In Fig. 21 we see that the pole approximation in Eg.( P03| ) which yields the 
power law Eg. ( ^04] ) is rather close to the exact numerical result in the full 
range of 5. The approximation Eg. ( |306| ) for the moment c*-^-* approaches the 
correct power for sufficiently small 5(ln'(9/5 > 1), but rises stronger towards 
larger 5 ^ 'd where the integral over fli approaches the singularity of F*^^^ at 
{}i = 0. This demonstrates the importance of the moments for the test 
of power behaviour. 

It may be surprising at first sight, that a power law is found for the cu- 
mulant moments and not for the factorial moments which are the integrals 
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of the density p'-"^ in Eq.( |179| ), while on the other hand p'^'^\'di2) is power 
behaved. In the Appendix D we perform some consistency checks proving 
that indeed this is the case. In particular we integrate successively the cu- 
mulant correlation r*^^)(i7i, ^2) up to the multiplicity moment /2 with either 

P^'^\'^i2) or h^'^\VL,5) and S'^\'d^5) as intermediate steps. 
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Figure 1: Diagrammatic representation of a) the integral equation (|T5|) for 
partonic densities and b) of its solution, Eq.(|6lD. 
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Figure 2: Kinematics of a QCD jet in the transverse plane. Allowed phase 
space of the intermediate parent K is shown as the shaded regions. The 
large circles indicate the cone of half opening B around the primary parton 
of momentum P; left column: (a) precise bounds for p'p in Eq.(|^), (b) for 
in Eq.([l|) or A^,^^ in Eq.(|lD, (c) for in Eq.®; middle and right 
column: integration areas in the pole approximation. 
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Figure 3: Kinematics of the inclusive production of a single parton. 



Figure 4: Second iteration of the integral equation ([TsD . The vertex Ki K2 
can be interpreted as the first term of the power expansion of the resolvent 




Figure 6: Definition of tlie relative angle 1^12 for particles in a forward cone 
around a primary parton of momentum P. 




Figure 7: Normalized 2-parton angular correlation r{'di2) in the forward 
cone of half opening angle = 1 for asymptotic and finite energies {Y = 
ln(P/go) = 5,go/A = 2). 




Figure 8: Normalized 2-parton correlation functions at finite energies for 
quark and gluon jets and in the infinite energy limit, (a) Differential nor- 
malization (|123|) and (b) global normalization ( |124| ). 
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Figure 9: Singularity structure of the terms building up the connected cor- 
relation function F*^"^ (for n=3). A line between particles i and j denotes a 
singularity 




Figure 10: Definition of angles in a sidewise ring and a sidewise cone refered 
to as 1-dim and 2-dim configurations respectively. 
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Figure 11: Maximal sidewise cone used to define h{6) in Eq. (|193|) . Kinematics 
corresponds to the subprocess (B) discussed in the Appendix J, i.e. P — > 
K -^3^ 1(c). 
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Figure 12: Normalised moments for (a) the sidewise ring (D=l) and (b) 
the sidewise cone (D=2) from Eg. ( pi6| ) using the high energy approximation 
(evaluated for momentum P = 45 GeV, = 30°). For e smaller than indi- 
cated by the arrows the asymptotic expressions apply, (c) Energy dependence 
of the normalized moments. 
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Figure 13: Rescaled 2-particle correlation f{s) vs. scaling variable e for dif- 
ferent primary energies P from the parton MC. The full curve represents 
the high energy limit of the DLA, the dashed curve the prediction for quark 
jets at P = 45 GeV; the normalization of the curves is adjusted. The insert 
shows the energy dependence of the same quantity at fixed e. 
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Figure 14: Rescaled two particle correlation as in Fig. 13 but with differential 
normalization as in Eq. (|123|) . The MC data don't approach the scaling limit 
represented by the curve. 




Figure 15: Cumulant and factorial moments C^^^ and F^^^ for different dimen- 
sions D and primary momenta P from the parton MC. The curves represent 
our results for C^^^ whereby the normalization is adjusted to fit the MC data. 
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Figure 16: Rescaled cumulant moments C2 as defined in Eg. ( |222| ) for tlie 
sidewise ring (ID) and tlie sidewise cone (2D) for different jet momenta P 
as function of e in comparison witli tlie fiigli energy DLA prediction 2/5(1 — 
uj{e,n)/n), normalized to zero at e = 0. 




Figure 18: Dependence on the order n of the D = 1 rescaled facto- 
rial moments. The curves represent the DLA predictions, normahzed to 

^^(O) = F('^)(0) = 0. 
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Figure 19: Comparison of results on correlations and moments obtained at 
the parton level and after hadronization for u-quark jets. Results for a su- 
perposition of all jet flavours produced in e'^e~ — > hadrons don't differ ap- 
preciably. 
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Figure 20: Angular distribution dn/drj for different jet momenta P with 
respect to the sphericity axis and with respect to the particle's directions 
weighted by their energies. The full (dashed) curve represents the exact 
(approximate) DLA results where the normalization is adjusted. 
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Figure 21: Numerical check of various approximations in moment calcula- 
tions: the straight line represents the function 7T6'^h{Q, 6) using pole approx- 
imation Eq.( p04| ) and the dots the result of exact numerical integration in 
Eq. (|303|) . The dashed curve represents the exact calculation of the moment 
c(2)(t9,5) in Eq.(|T80l), for Y = \nP/Qo = 5 and 6 = 1. 



